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DECISION MAKING AND NEW MATHEMATICS! 


Gerald L. Thompson* 
Dartmouth College 


Anyone familiar with the history of science is aware that the development of physics 
and the development of mathematics are intimately connected. For example, Newton's 
Principia Mathematica contained, besides the physical theories with which his name is most 
often connected, also the origins of calculus and some very interesting investigations of 





algebraic curves (pure mathematics). Many other examples of this influence can be cited. 

It also is the case, although less well known, that problems of decision making have changed 
the course of mathematics. I am not prepared to give a historical study of this influence but 
rather to discuss some current problems of decision making, stimulated by the social or 
behavioral sciences, in particular economics, as they are influencing the creation of new 
mathematics today. 

What I want to do here is to discuss a dozen or so new kinds of problems, arising 
directly or indirectly from decision making questions, which have been studied by one or more 
mathematicians. I will accept the assertion that if some new mathematical theorems have been 
discovered as the result of such studies, then these new problems are changing the course of 
mathematics. 

Surprisingly enough, the problems to be discussed are of very recent origin. Indeed, 
the last fifteen years suffice to cover their publication dates, although their roots in many 
cases go back much further. 

It is not possible for me to go into the technical details of the problems here. But I can 
illustrate them by means of simple examples and in this manner present an idea of the wide 
range of decision making problems that are being studied at the present time. 

The first class of problems to be discussed are the so-called ''games against Nature." 
A game against Nature is basically a statistical problem which I illustrate by means of the 
following oversimplified example: Suppose that you are a farmer and must choose in a certain 
year whether to plant wheat or corn (you can plant only one of these). Suppose also that hot, 
rainy weather is good for corn and bad for wheat, and that cool, dry weather has the reverse 
effect (and these are the only possible kinds of weather). To be specific, let's assign some 
humbers to illustrate these effects, as in the table below. 


*Manuscript received 5 October 1956 


4 lPaper was given at the Conference on Methods in Philosophy and the Sciences in New York 
City, May, 1956. 
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Weather (Nature) 
Hot, Rainy Cool, Dry 


Corn 10 5 
Wheat 3 8 




















The numbers represent (say) the yield in dollars per acre under the conditions specified. 
Observe that corn is, in general, a more profitable crop than wheat but that, in a bad year for 
corn, wheat is more profitable than corn. How shall you (the farmer) act, assuming that you 
have absolutely no information as to how Nature will choose the weather in the coming year 
(there is no long range weather forecast, etc.). 

There are essentially two different kinds of solutions that have been proposed by 
various people for the problem: (a) Make an assumption as to what Nature will do and then 
maximize your own position relative to this assumption. (b) Try to choose your course of 
action (strategy) so that, regardless of what Nature does, you have in some sense protected 
yourself against the worst that could happen. 

As representative of the decisions of type (a) we mention the Laplace decision rule: 
Choose whichever row has the largest average. A probability rationale for the Laplace rule 
can be given as follows: in the absence of any information as to what probability Nature has 
placed on the kinds of weather, assume that they are equally likely; e.g., in our example above 
assume that the probability of hot, rainy or cool, dry weather is each equal to one half; then 
compute the mathematical expectation* that results from each course of action under this 
assumption and act to maximize this expectation. 

In the example above, the farmer's mathematical expectation for a given course of 
action, that is, for the choice of a given row, is the average of that row. The first row (plant 
corn) has the larger average, namely, 7.5; hence our decision, using the Laplace criterion, is 
always to plant corn. 

The Wald and Savage decision rules are of type (b). Let us consider the Wald or 
minimax rule first. It is the following: Assume that the table of payoffs is that of a two-person 
zero-sum game and use as a decision rule an optimal strategy for that game. In our example 
above, the optimal strategy that the farmer should use is to choose each row with probability 
one half; for example, he could flip a coin and choose the first row if (say) heads turns up and 
choose the second if tails turns up. If the farmer uses this strategy year after year, then 
about half the time he will plant corn and half the time wheat. In this case, regardless of what 
weather Nature actually chooses, the farmer's expectation is 6.5, i.e., his average earnings 
over the years are 6.5 dollars per acre. 

The Savage or minimax regret rule can be described as follows: Suppose that you knew 
what Nature's choice of a column would be; then define as your regret for having chosen a given 
course of action, that is a given row, the difference between the entry you receive in the indi- 
cated column and the maximum entry in that column. For the example above, the table of 
regrets is 




















*To explain the concept of mathematical expectation which we shall use frequently in this paper, 
suppose that there are n events of which one, and only one, can happen; suppose they give to us 
payoffs (monetary) inthe amounts aj, a2, ..., an; suppose that we know there are probabilities 
assigned to the occurrence of these events, say, Pp], p2,...+» Pn» Where p; 2 0 and p) + p2 + 
+. .+ py = 1; then our mathematical expectation E is 


E = a)p, + a,p,+...+ app 


n° 
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The minimax regret rule now is: Regard the regret table as a two-person zero-sum game 
and use as a decision rule an optimal strategy in that game. For our example the minimax 
regret rule gives the result: Plant corn with probability 0.7 and wheat with probability 0.3. 

Observe that each of these decision rules gives a different answer. There are several 
other decision rules which give still other different answers. Each of these rules satisfies 
certain desirable properties, but none satisfies all desirable properties. For further discus- 
sion we refer the reader to Milnor's paper [9], on which the present discussion has been based. 

The Wald decision rule which we have illustrated for a finite matrix game has much 
more general interpretations in the Theory of Statistical Decisions. There a pure strategy of 
the experimenter is to choose a given decision rule as a method of testing a hypothesis, and a 
pure strategy of Nature is the choice of a distribution. Relative to a choice of a decision rule 
by the experimenter and a distribution function by Nature there is a payoff, called the risk, 
that gives the cost of experimentation for that case. Then the decision problem is reduced to 
an infinite matrix game where Nature has at most a countable number of pure strategies. An 
optimal strategy in this game provides Wald's solution to the statistical decision problem. 

Continuing our discussion, let us briefly consider the dynamic programming problems 
of R. Bellman [2]. Here an example will suffice to state the problem. Suppose that we are the 
owner of two gold mines of which one is a surface mine yielding very low grade ore and the 
other is a deep mine yielding high grade ore. The surface mine is easy and inexpensive to 
work but has low yield, while the deep mine is difficult and expensive to work but has high 
yield. Suppose that we have a fixed amount of capital which we desire to spend on the mining 
operations, At the beginning of each time period, which may be a month or a year, etc., we 
divide the available capital into two parts and use one part for surface mining and the other 
for deep mining. At the end of each time period a certain part of the capital remains for use 
during the next time period. The problem is to decide how to allot the portions of capital so 
that our returns are maximized. 

The main difference between this and the preceding problems is that the present prob- 
lem involves decisions made repeatedly over several time periods. Possibly even an infinite 
sequence of decisions is needed, each depending in some way on past results and an estimate 
of what will happen in the future. The most general such problem that could be stated is 
insoluble, but Bellman has shown the existence of optimal decision rules under certain sim- 
plifying assumptions. 

The problems discussed so far belong to a class called one-person games. Thus in 
each of them there was only one decision-maker who controlled some, but not all, of the vari- 
ables in the problem. The rest of the variables were controlled by Nature, about whose actions 
usually no a priori assertion could be made. Thus Nature could or could not be a malevolent 
participant in the game. The most conservative or pessimistic assumption about her behavior 
was that she was malevolent, and that assumption lead to minimax decision procedures. More 
Optimistic assertions about Nature's behavior lead to other, less pessimistic, decision rules. 

We now turn to the class of so-called two-person zero-sum games in which there are 
two active participants. The initiators of this, as well as of most other parts of the theory of 
games, are von Neumann and Morgenstern [10]. The minimax theorem of the theory of games 
assures that a solution, in the form of optimal (mixed) strategies for each player and an 
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(expected) value for the game, exist in all cases. Since there is some confusion as to the 
meaning of the minimax theorem, let me illustrate it by means of the customary example of 
matching pennies. The diagram below, where we have named the players Jones and Smith, 
gives the explanation. 


Smith 
Head Tail 


Head 1 -1 
Tail -1 1 





Jones 

















Here the entries in the table are the payoffs which Jones receives from Smith (we assume that 
Jones is trying to match Smith). The value of the game is zero, and optimal strategies for 
either player are to choose head or tail with equal probability, i.e., with probability 1/2. How 
can we interpret the value (zero) of the game? Observe that in a single play of the game each 
player will either win or lose one penny, whereas the value of zero might seem to indicate that 
neither player could either lose or win. Here we must remember that the word value means 
expected value in the sense of probability. Thus if the game is played many times one can say 
that with high probability the average earnings of each player become close to zero. In a single 
play of the game one can say that the optimal strategies make the probability of loss equal to 
the probability of win. 

It is well known that players of this game realize their optimal strategy by "flipping" 
their coin, that is, leaving their choice of head or tails to a chance device that chooses head 
and tail with probability 1/2. Thus the player himself does not know beforehand what choice 
he will make. Here is a case where spying would help. For example, if Jones could find out 
what choice Smith is going to make he could adjust his own choice accordingly and always win. 
Smith can prevent Jones from finding out (short of some psycho-kinetic effect) by leaving his 
choice up to the chance device and, in fact, not knowing himself exactly what choice he will 
make. 

The idea that "playing a game optimally" should involve "leaving your choice up to a 
chance device" seems flippant to many people, particularly when game theory is applied to 
more important kinds of games than matching pennies. Another common criticism of the 
results of the minimax theorem is the following: Suppose that Smith is stupid enough always to 
choose the same side of his coin, say heads. In the course of playing the game several times 
Jones would certainly observe this behavior and, one would think, should take advantage of it by 
always himself choosing heads. But the latter behavior is not optimal! Indeed, if Jones plays 
optimally, then his expectation remains zero regardless of what Smith does, and in effect he is 
protecting Smith from the follies of his own (Smith's) stupidity! 

The answer to the first criticism is that the minimax theorem proves that no other 
behavior can do better in all situations than the optimal behavior. It shows that the flipping of 
a coin, or the use of similar chance devices, is actually a serious and important way of making 
strategic decisions. The answer to the second criticism can be made in two ways. Some 
theorists would state that in order to take advantage of an opponent's mistakes it is necessary 
to deviate from an optimal strategy and take risks, and risks of that sort are not a proper 
subject of study in game theory. Another answer might be that if Smith plays in the stupid 
manner described above and Jones finds out about it, then Jones can suitably redefine the game 
for himself (e.g., cross out the second column) and play minimax in the redefined game. The 
optimal strategy in the new game is the one suggested above. 
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The theory of linear programming was developed by G. Dantzig and others about 1948. 
The classical example of a linear programming problem is the diet problem. Everyone must 
eat to live and must obtain minimal requirements of certain nutrients to remain healthy. Also, 
most people have limited incomes and hence wish to satisfy their diet requirements while at 
the same time minimizing their grocery bills. The amounts of nutrients contained in foods 
varies linearly with the quantity purchases and, if the quantities are small, so does the cost of 
the foods. Thus the problem is: What should our diet be in order that we obtain our minimum 
diet requirements and pay least for groceries? If the palatability of our food also is important 
we can add the requirement that at least certain minimum amounts of certain foods be eaten. 

Many other examples of linear programming problems can be cited from applications 
to military, economic, and business problems. 

The solution of a linear programming problem can be shown to be essentially equivalent 
to the solution of a two-person zero-sum game. The mathematical definition of a linear pro- 
gramming problem is that it requires the maximization or minimization of a linear function of 
n variables whose domainis a convex set of values. Because solutions of such problems occur 
on the boundaries of the convex set, it is not possible to use the ordinary techniques of calculus 
to find them. Instead, finite methods that can be programmed on high-speed electronic com- 
puters are used. 

As another application of two-person zero-sum game theory to decision problems, I 
would like to mention some very recent work by John Kemeny, Oskar Morgenstern, and 
myself [7]. We considered a generalization of the von Neumann model of an expanding economy 
and found that it was intimately connected with the theory of games. It turned out that, although 
the building blocks of the model were ordinary economic concepts such as inputs, outputs, 
interest rates, expansion factors, etc., the type of mathematics needed to analyze the model 
and solve the problems connected with it was the theory of games. In other words, although 
game theory was not used at all in the formulation of the model, the mathematics that had been 
developed for game theory turned out to be precisely the kind of tool needed for its analysis. 

As far as I know this is the first time that the mathematics of game theory has been so applied. 
We found that in the generalized von Neumann model, which is a very general type of economic 
input-output model, there can be at most a finite number of expansion rates at which the econ- 

omy can expand in equilibrium, and for each of these expansion rates, the interest rate is equal 
to the expansion rate. 

We turn now from the theory of zero-sum two-person games to nonzero-sum two-person 
games. In so doing we leave a well-ordered subject that is dominated by an intuitively accept- 
able theorem, the minimax theorem, and approach a subject in which there is no widely accept- 
able or accepted theory. In fact, the subject is such that there appears to be little hope at 
present of finding a universally acceptable theory. 

The only general theorem that has been proved in this area is the one due to J. Nash 
concerning the existence of equilibrium points. Briefly, it says that in any two-person game 
there always exists an equilibrium point, that is, a pair of strategies, one for each player, 
having the property that neither player can improve his position in the game by changing his 
own strategy while the other player keeps his strategy fixed. In other words, at an equilibrium 
point neither player can by his own efforts improve his position. Actually the idea of an 
€quilibrium-point solution goes back quite far in the economic literature—at least as far back 
as Cournot. There are, however, objections to equilibrium points as solutions of games, as 
the next example will illustrate. 
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Suppose that there are two entrepreneurs, specifically let us say gasoline station oper- 
ators, whose places of business are located on opposite sides of the same street. We may as 
well call them Jones and Smith, following our previous example. Suppose that all other 
gasoline stations are located at least several blocks away from Smith and Jones so that their 
only immediate competition is with each other. We shall assume that there are only two prices 
which each can charge for his gasoline, call them high and low. At the beginning oi 2ach day 
they must decide which price they will charge for the day and they cannot change price during 
the course of the day. The strategic decision each must make then is the price he will charge, 
Let the results of their decision be given in the diagram below: 


Smith 
High Low 
High (10,10) (6,16) 
Low (16,6) (7,7) 








Jones 














In each pair of numbers enclosed in the parenthesis we indicate by the first entry the amount 
that Jones gets and by the second the amount that Smith gets. Let us see the rationale behind 
these numbers. If both charge the high price then each receives 10. If both charge the low 
price, which we can interpret as a "gas war," then, even though the volume of business may be 
greater, the low price cuts their income to 7 for the day. Finally, if one charges the high price 
and the other the low, the low-priced man can draw business not only away from his immediate 
competitor but also from elsewhere; and the volume of his business yields him 16 while the 
high-priced man gets only 6. Whether or not you agree that these numbers are realistic, I 
would like to keep them as chosen to illustrate several points about nonzero-sum games. 

The equilibrium point in the game is (7,7), as can easily be checked. At this point each 
player charges the low price, i.e., the gas war prevails. If either player changes his strategy 
to that of charging the high price, he can only decrease his own income. To most people the 
solution seems irrational, since at the (10,10) point each player is better off than at the equi- 
librium point. But the (10,10) point is not stable, since either player can, by switching to the 
low price, improve his take from 10 to 16 (assuming that his opponent does not also change). 
It is also true that at the (6,16) or the (16,6) point the total amount they take in is greater than 
the total amount they take in at the (10,10) point, and so either of these points is preferable to 
the (10,10) point. Yet, again, neither of these two points is stable in the above sense. 

The solutions that have been proposed for the game are the following: We have already 
mentioned the equilibrium point (7,7). The following cooperative solution has been proposed: 
Let one player, say Jones, always charge the high price and let Smith always charge the low 
price. Then Jones gets 6 and Smith gets 16 and their total take is 22. Now equalize the situa- 
tion by having Smith make a side-payment of 5 to Jones and each gets a total take of 11. If 
there is some reason, moral or legal, why such a side-payment should not be made, then the 
players could adopt the following strategy: On alternate days let Jones charge high and Smith 


low and on the other days let the reverse be true. Then the average payment to each will be 
11. 


The scheme just described for mediating or arbitrating this two-person game is due to 
Nash and Raiffa and is also one of cooperative solutions proposed in the von Neumann- 
Morgenstern theory. It has been used in a very charming way by an English philosopher 
Braithwaite [3] to resolve another, somewhat similar situation, which he chooses to call a 
moral situation. 





DECISION MAKING AND NEW MATHEMATICS 147 


There is still another solution to this game which comes out of some work being done 
by M. Shubik and myself. Suppose that Jones is richer than Smith and suppose that at the 
equilibrium point (7,7) both players actually lose money. Since this game is played over and 
over again every day, we could imagine Jones making a threat to Smith as follows: "If you do 
not join with me in charging the high price, the day after you lower your price I shall lower my 
price and keep it down until you are bankrupt."" Since Jones has more capital than Smith, the 
threat is effective and, if carried out, would result in the closing of Smith's business. Suppose, 
on the other hand, Smith tried to make a similar threat. Since he would go out of business first 
if he carried out the threat, it is not effective for him but merely suicidal. In this case Jones 
is able to enforce the (10,10) point and has the dominant position in what otherwise appears to 
be a symmetric game simply by being richer. We are continuing our work on the problem 
along this and other lines. 

There are cases of three-and-more-person games which give illustrations as to the kinds 
of results that are being obtained there, but I will leave the theoretical realm and enter into 
the experimental realm by discussing some very unusual experimental results recently 
obtained by psychologists. These results will have a direct bearing on the concept of ration- 
ality, and I shall have a few comments about that concept later. 

Let us digress a moment and talk about baseball. Suppose that, in a certain year, the 
following situation arose: The Yankees beat the Dodgers, the Dodgers beat the Giants, and the 
Giants beat the Yankees. Most probably the man in the street or the sports commentator 
would say that this was an "upset.'' It is also an example of an intransitive triple, namely, 
given three events, say A, B, and C, and an ordering relation ''greater than" if it happens that 


A is greater than B, B is greater than C, and C is greater than A, then we say that A, B, and 
C form an intransitive triple. Similar intransitive relations also occur in the "pecking order" 


of chickens in a barnyard. 

Most armchair thinkers about rationality assume that a "rational" person is consistent 
in the sense that his preference orderings for goods and services are transitive. Recent 
experimental work by Flood [6] and Papendreou and theoretical work by Arrow [1] and May [8] 
have shown that when groups of people are concerned, there is a very real possibility of an 
intransitive ordering being set up by them. Even a single individual may have intransitive 
preferences, and I think that economists will soon have to include this possibility in their 
theories, 

Another place where apparently irrational behavior on the part of individuals and 
groups of individuals has been observed is in the area of psychological learning theory. I shall 
describe an experiment from this area, one which I think has the most intriguing results of 
recent experimental psychology. Suppose that the subject (who may be a rat, a goldfish, a 
pigeon, or a human being) is placed in a room in which there are two "slot machines" and is 
instructed to play them. The experimenter has adjusted the machines so that one of them pays 
off 50 percent of the time and the other pays off never. What will the subject do? The experi- 
ment has been performed many times in many different guises and the results are essentially 
the same, Asymptotically it is observed that the subjects will tend to choose about 2/3 of the 
time the machine that has the 50-percent payoff and about 1/3 of the time the machine that 
never pays off. This result seems quite paradoxical since the game theory solution to the 
problem is always to choose the machine with the higher probability of payoff (after a few 
observations have established which is higher). Here is a case in which experiment gives 
different :esults than pure thinking. Learning theories which predict the observed result have 
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been devised by Bush and Mosteller [4] and Estes [5]. The mathematical study of these learn- 
ing models has given rise to difficult mathematical problems, some of which are presently 
unsolved. 

I can state another, to me astounding, result from experimental psychology, after a 
short digression. A man in the street, as well as the German philosopher K. Marbe, believes 
in the so-called gambler's fallacy, namely, that if in the repeated tossings of a true coin a 
long run (say 3 or 30 or 30,000,000) of heads have turned up, then the probability of a tail 
turning up the next time is larger than 1/2. Experimental evidence does not confirm this 
belief. But if one takes a subject who indicates that he consciously believes in the gambler's 
fallacy and places him in an experimental situation of the kind described above, then, accord- 
ing to Estes, he will produce, in the course of his experimental run, a very nice random 
sequence—one containing no hint of the gambler's fallacy. We thus have the paradox of a 
subject who consciously disbelieves in random sequences, but who can produce, while appar- 
ently doing something else, very good random sequences. 

I make the last two observations because I believe that in economics, and in other 
decision making fields, the day of armchair thinking about the rationality of the decision- 
maker is at an end. In attending meetings of professional economic societies I have observed 
that there is some reluctance on the part of economists to accept and include in their theories 
such empirical findings as the ones mentioned above. In my opinion, such results cannot be 
ignored and will have important influence on future economics. 

In conclusion, let me say that in the course of this rapid survey of decision making I 
have mentioned several theories and indicated the kinds of mathematical problems to which 
they lead. In most cases, entirely new kinds of mathematics had to be developed and were 
developed by the cooperation of mathematicians and social scientists. There is no reason to 
think that this is a temporary trend or fad. On the contrary, I believe that it will continue to 
be of importance in the future. 
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A PLAN TO ALLOCATE AND PROCURE ELECTRONIC SETS BY THE USE OF 
LINEAR PROGRAMMING TECHNIQUES AND ANALYTICAL METHODS 
OF ASSIGNING VALUES TO QUALITATIVE FACTORS 


Jack W. Smith* 
Navy Management Office 
Department of the Navy 


I, INTRODUCTION 

The Bureau of Ships must allocate limited available quantities of electronic sets to 
ships of the fleet. Accordingly, large-scale allocation plans are developed to specify which 
ships will receive new sets and which ships will not. The development of such allocation plans 
requires simultaneous consideration of many factors. Some of these factors are: 

(1) The priority of ships and requirements; 

(2) The effectiveness of sets already installed; and 

(3) The importance of taking full advantage of all of the features built into each of the 
new sets. 

The problem of allocating is further complicated by the fact that often a particular kind 
of set can be used to fill many different requirements while a particular requirement can be 
filled by numerous different kinds of sets. Furthermore, the interchangeability of available 
sets sometimes makes it necessary to consider simultaneously as many as 1000 different 
Ships, having up to 50 different kinds of sets installed, and 10 different kinds of sets available 
for allocation in quantities ranging from 1 to 300. 

Until recently, intuitive methods have been used exclusively in developing these plans; 
there was no other way to determine if available sets were being utilized in the best possible 
way. Now a method has been devised in which linear programming techniques are used to 
develop allocation plans, and the method is being extended into the field of procurement 
Planning. In addition, analytical methods have been found for assigning quantitative values to 
factors or concepts which are basically qualitative in nature. | 





x 
Manuscript received 15 January 1956 


For their help in the evolution of this method, the writer is indebted to Woodbury, 
Cannon, Heller, McShane, and Stephenson of the George Washington University Logistics 
Research Project; Polachek and Suzuki of the David Taylor Model Basin Applied Math 
Laboratory; Weyl of the Office of Naval Research; and Kruskal and Aumann of the 
Princeton Analytical Research Group. 
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Il. BACKGROUND INFORMATION 

The Bureau of Ships procures and allocates electronic sets to fill ships’ requirements 
specified by the Chief of Naval Operations (CNO). . These requirements are essentially func- 
tional descriptions of what an electronic set should do. For instance, CNO might state that 
each of a certain group of ships requires radio transmitters covering a frequency range 
from 20 to 30 Mc, having a reliable range of 300 miles when received by an average shipboard 
receiver, and capable of both continuous-wave and amplitude- modulation emission. 

A model is a group of essentially identical sets built under the same specification. In 
general, any model which will meet the criteria specified in the requirement and which is 
within certain limitations as to weight and size is acceptable. There may be a number of such 
models in stock or under contract. Usually one of these is considered to be the preferred 
model for a specific requirement on a specific ship. The capability of a model may overlap 
different requirements. For instance, a 10- to 30-Mc transmitter would fill both a 10- to 
20-Mc requirement and a 20- to 30-Mc requirement. This same transmitter would be a 
partial substitute for a 20- to 35-Mc requirement. 

As used in this paper, a ship class is a group of ships having identical requirements 
and missions. The individual ships within a class may be quite dissimilar in that their require- 
ments have been filled with different models. Also, different requirements may have been 
filled or partially filled on some of the ships within a class and not on others. 

An allowance is the number of times a specific requirement should be filled on each 
ship within a given ship class. This allowance may be as great as 15. In such cases of "mul- 
tiple allowances" it is usually considered more important to give a ship its first set than it is 
to give another ship in the same class its fifteenth, or even second, set. 


Ill. PRINCIPAL FACTORS INFLUENCING ALLOCATION DECISIONS 
The Priority of Requirements and Ship Classes (The Material Improvement Plan) 

Because of the limitation of procurement funds and changing requirements, the total 
available number of preferred sets and suitable substitutes combined often is not sufficient to 
outfit all ships. Furthermore, installation funds often are limited. To guide the materiel 
Bureaus in the allocation of limited sets and limited installation funds, CNO provides the 
Materiel Improvement Plan (MIP). This document lists groups of requirement-ship class 
combinations in priority order. Within each requirement-ship class group, ship classes are 
listed in order of importance. A requirement such as surface search radar might be group 
priority number one for certain ship classes and group priority number sixteen for other ship 
classes. Thus each requirement-ship class combination has its own unique place in the over- 
all priority structure indicated by the MIP. 

The form of the Material Improvement Plan is as follows: 














Priority Requirement Ship Class 
1 Surface Search Radar CVA 
DDR 
CA, etc. 
VHF Transmitters DDR 
DM 
APD, etc. 
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Priority Requirement Ship Class 
3 Portable Transmitters DDE 


DE 
DER, etc. 


Surface Search Radar AGC 
APA 
BB, etc. 


The MIP does not take into consideration the relative importance of completing a multiple 
allowance which is partially filled, as opposed to supplying a ship its first set. Nor does the 
MIP consider the degree of effectiveness of sets already installed. 


The Effectiveness of Sets Already Installed 
To fill a particular requirement, many different models may have been installed. These 
installed models may be perfectly satisfactory, absolutely unsatisfactory, or fall anywhere 
between these extremes. The degree of effectiveness depends upon the following: 
(1) The Requirement 
A model may not contain all required features even when it is operating at peak 
performance. 
(2) Anticipated Maintenance Difficulty 
With different models different degrees of maintenance difficulty (down time, etc.) 
are to be expected. 
The Ship Class 
The preferred model is often the heaviest one which can be installed on a particular 
ship class. Thus the small surface search radar preferred for a small ship may be 
considered substandard for a larger ship which can carry a heavier, more powerful 
set. 


The Effectiveness of Sets Available for Allocation 

There may be a number of different models available which could be applied to a par- 
ticular requirement. Furthermore, there may be a number of different requirements to which 
an available model could be assigned. A typical problem might involve a structure of require- 
ments and available models as follows: 





Ship Class a ' Ship Class b 
Available Models, in , Available Models, in 
Order of Preference aepaeias Order of Preference 
1 A,C,B 1 B, A, C 
C,D,E 2 D,C,E 

E E 





Requirement 























Requirements 1 and 3 have no applicable available models in common. Requirement 2, how- 
ever, has applicable available models which are also applicable to requirements 1 and 3. 
Accordingly, all three requirements and all five models must somehow be considered simul- 
taneously when an allocation plan is developed. 
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As in the case of installed sets, the applicability of sets available for allocation 
depends upon the requirement, the anticipated maintenance difficulties, and the ship class, 
The way the ship class becomes a pertinent factor is illustrated in the above example by 
showing a switch in the order of preference of models as the ship class changes. Model A is 
the preferred model for requirement 1 on ship class a. Model B is the preferred model for 
requirement 1 on ship class b. This could happen, for instance, if model A was quite large 


and capable of the greatest range and ship class a was not weight or space critical while ship 
class b was critical. 


Utilization of Features 

Consider models C and E in the typical requirement-available model structure. Each 
of these models is applicable to two different requirements. This can occur in a number of 
ways. Requirement 2 might simply be requirement 1 plus some additional feature such as 
increased range, automatic self-adjustment, etc. If model C had all features included in 
requirement 2, it would automatically be capable of filling requirement 1, providing additional 
weight and space were not overriding considerations. But if model C were allocated to fill 


requirement 1, some of its features (i.e., those which are necessary only for requirement 2) 
would be wasted. 





With model E the situation would be different. Suppose that some features necessary 
to fill requirements 2 and 3 overlap.Model E might fill requirement 3 perfectly, with none of 
its features wasted. Model E would then be a substandard substitute in requirement 2 and, in 
addition, would waste some of its features in this application. 


Time Phasing 

Different models are produced and delivered at varying rates. Quantities of some 
models applicable to a particular problem might already be in stock. Others might be on order 
and scheduled for delivery. The further into the future delivery is scheduled, the less accu- 
rately can the date of delivery be predicted. Generally, deliveries estimated at 18 months in 
the future can be expected anywhere between 12 and 24 months. 

"Firm" ships’ overhaul schedules are received annually. Each schedule covers an 
entire year and is received approximately 4 months before the first ship in the schedule is to 
commence overhaul. The schedule changes somewhat during the year. Most changes involve 
advancing or delaying the overhaul of a ship one to three months, although ships are occasion- 
ally added or deleted. A "tentative" overhaul schedule covering the year immediately following 
that covered by the "firm" schedule is received with the "firm" schedule. 

Consideration of time in the allocation of electronic sets leads to some interesting 
questions. For instance, 

(1) Should a set be held for a high priority ship even if the ship's predicted availability 

is quite far in the future? In other words, is time phasing necessary at all? 

(2) Is time phasing for more than one year in the future worthwhile, in view of the 

relatively poor long-range predictability of equipment delivery schedules and ships' 
overhaul schedules ? 


IV. ALLOCATION PLANS PREPARED BY FORMER INTUITIVE-JUDGMENT METHODS 
CNO requires allocation plans to be prepared once a year for all ships scheduled for 
overhaul during a future two-year period. The plans become the basis for authorizing 
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installations, planning and scheduling shipyard work, distributing sets in stock or delivering 


from contract, and preparing budget requests for installation funds. 

The cost of electronic sets allocated in these plans is often more than $100,000,000. 
Records of sets installed are maintained in a deck of 500,000 IBM cards. The allowances of 
all ship classes are furnished by CNO. Production schedule estimates, stock records, and 
overhaul schedules are maintained. Using this information, a group of Supply Requirements 
Specialists and Electronics Engineers prepares allocation plans. 

When there is a shortage of available sets and there are no other competing factors, 
the Material Improvement Plan priority structure controls allocation. Available sets are 
allocated down the priority list as far as they will go. 

When there is a shortage of sets and multiple allowances are a pertinent factor, interim 
(or reduced) allowances are sometimes developed to enable assignment of at least one set each 
to some of the lower priority ship class-requirements. Then the available sets are allocated 
as far down the reduced-allowance priority list as they will go. 

When there is a shortage of sets, and the sets are applicable to more than one require- 
ment, and there are installed sets which fill requirements in varying degrees of effectiveness, 
and there are multiple allowances on some ship classes, etc., solutions are developed by the 
use of intuitive methods and simplifying assumptions. Sometimes installed models of border- 
line effectiveness are assumed to be either completely satisfactory or completely unsatisfac- 
tory. If an available model can fill two requirements, the decision may be to use it in only one 
of its possible applications. It is sometimes decided to let the effectiveness of installed sets 
be the primary and controlling consideration while the MIP priorities are the secondary con- 
sideration, or vice versa. Other factors are handled by making other simplifying assumptions 
and by developing procedures for the relation of the factors. For instance, it might be decided 
to use available model A sets to replace all installed model F sets in requirement R, and then 
to use the remaining model A sets to replace installed model G sets in requirement Ro, in each 
case following the MIP priority structure. 

The writer, who was responsible for supervising the development of electronics allo- 
cation plans, never felt confident of the results. There were often too many ways that the 
simplifying assumptions could have been made. Also, there were too many procedural methods 
possible, and each would lead to a different allocation plan. 





V. FORMULATION OF THE PROBLEM 

In view of the significance of the allocation function, the abundance of basic data, and 
the questionable results of intuitive methods, it was decided to investigate the possibility of 
developing analytical methods for preparing allocation plans. To this end a "model" of the 
problem was developed. The basic idea of this "model" was that values or "weights" could be 
assigned to various factors and that by combining the weights a numerical value could be found 
for each possible way to assign each available equipment model. Thereafter, the method of 
assigning equipment models could be varied until the sum of weights corresponding to the 
assignments was maximized. 

The procedure was to be as follows: 

Step (a). Prepare tables for each ship class-requirement in which all equipment 
models, both installed or available for allocation, are arranged in preference order. 

Step (b). For each preference table assign "Improvement" weights corresponding to the 
improvement gained by replacing any equipment model with each model higher in the preference 
table. For instance, if in the preference table for a particular ship class-requirement model A 
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was preferred over model B, which in turn was preferred over model C, then three different 
improvement weights would be assigned: one each for A replacing B, A replacing C, and B 
replacing C. The method by which these Improvement Weights are determined is explained in 
Section VI. 

Step (c). Assign MIP Weights corresponding to the priority of each ship class- 
requirement, as indicated in the Material Improvement Plan (MIP). 

Step (d). Considering specific ships and the specific models installed, combine the 
applicable Improvement Weights (step b) with the MIP Weights (step c) to obtain new weights 
which consider both factors. 

Step (e). Assign to the models available for allocation Utilization of Features Weights 
corresponding to the percentage of features utilized in each of their possible applications. For 
instance, a weight of 10 might indicate that none of a model's features are wasted, while a 
weight of 8 would indicate that approximately 0.2 of its features are wasted. 

Step (f). Combine these new Utilization Weights with the combined MIP-Improvement 
Weights (obtained by step(c)) to obtain a weighted list which contains each possible application 
of each model available for allocation. (One requirement-model installed combination on one 
ship would be listed a number of times and with different weights, since it would be possible to 
use a number of different models for replacement.) 

Step (g). Allocate available models (allowing one requirement-model installed com- 
bination on one ship to be filled only once) so that the sum of the weights is maximized. 

Step (h). Prepare a complete solution, ignoring the multiple allowance factor, i.e., 
considering the fifteenth equipment as important as the first. Thereafter, reduce all allowances 
by 24% (1 remains 1, 2 remains 2, 3 becomes 2, etc.) and obtain a new solution. If there is still 
a significant shortage, reduce all allowances by 48% (1 remains 1, 2 becomes 1, etc.) and obtain 
a third solution. Etc. 

Solutions based on ultimate allowances, and also on different kinds of interim or reduced 
allowances, would then be available. 

Step (i). Send all solutions or suitable abstracts to CNO for decision. (Steps (h) and (i) 
substitute a procedural-judgment method for a weighting system insofar as the multiple allow- 
ance factor is concerned.) 

Suzuki of the David Taylor Model Basin Applied Mathematics Laboratory and 
Stephenson of the George Washington University Logistics Research Project formulated the 
model as a transportation type linear programming problem, thus providing a systematic 
procedure for obtaining solutions which are optimal within the framework of the model. In this 
formulation the "worth" or "value" of assigning a specific equipment model to replace a specific 
installed equipment model in a particular ship class-requirement combination was considered 
as a function of three factors, namely, MIP Weight, Improvement Weight, and Utilization of 
Features Weight.2 These three concepts are basically qualitative, however whereas, for the 
transportation model, quantitative values are required. 

















2Time phasing of equipment delivery and overhaul schedules can be brought into the frame- 
work by expanding the matrix and at the sarne time discounting allocation worths corresponding 
to the time during which equipment must be held in storage. Current solutions, however, ignore 
the time phasing aspects since it is the normal procedure for equipment to be held in storage 
until the proper vessel is overhauled. The multiple allowance factor can still be considered by 
the procedural method described in steps (h) and (i). 
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vl. ASSIGNMENT OF VALUES OR WEIGHTS BY USING RANK ORDERING DECISIONS 
AND DEVELOPMENT OF THE FUNCTION BY WHICH MILITARY WORTH IS 
DETERMINED 

Generally speaking, it is possible to obtain agreement and authoritative decisions on 
rank orderings, i.e., A is better or more desirable than B, etc. The MIP priority structure 
provides such a rank ordering of requirements and ship classes. But it is difficult or even 
impossible directly to obtain agreement or authoritative decisions on quantitative values or 
weights, i.e., A is 2.3 times as important as B, etc. Even if such agreement could be obtained, 
the significance of the values would be questionable and it would be difficult to determine how 
to use them to arrive at an "over-all military worth" value. Accordingly, rank-ordering deci- 
sions have been used as a means of obtaining quantitative values or weights. 


The Assignment of Improvement Weights 
To arrive at Improvement Weights the writer first obtained technical agreement and 

authoritative approval of the rank order of the effectiveness of the various models available 
for or installed on a particular ship class for a particular requirement. This was called the 
"goodness" rank order and was of the form: 

Model A is more desirable than 

Model B, which is more desirable than 

Model C, which is more desirable than 


Model G, which is more desirable than 


Model H or a void (unfilled requirement). 
Next, agreement was obtained as to the rank order of the improvements gained by 
replacing any model with a better model, i.e., 





A replacing H provides a greater improvement than does 

B replacing H, which provides a greater improvement than does 

A replacing G, which provides a greater improvement than does 

B replacing G, which provides a greater improvement than does 

« g Gee. 

This improvement rank-ordering contains some very significant information. Obviously, 
the greatest improvement is obtained by replacing the worst model or a void by the best model 
(A replacing H in the example). But the next most desirable improvement could be either B 
replacing H or A replacing G. The decision depends on whether the difference between A and B 
or between G and H is greater.3 : 

After these rank orders were obtained, the writer found and assigned "goodness values" 
to each of the models in the goodness rank-order in such a way that they satisfied the following 
conditions: 

(1) The "goodness" values were consistent with the established "goodness" rank-order, 
i.€., goodness value for Model A > goodness value for Model B>... . > goodness value for 
Model H. 

a 


3The improvement rank-order agreed to in any specific case is only one of a very large 
number of different possible rank orders which are self-consistent and consistent with the good - 
ness rank-order as well. Such consistency seems to be inherent in reasonable improvement 
rank-orderings. 


419839 O - 57 -3 
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(2) The arithmetic differences of the goodness values were consistent with the estab- 
lished improvement rank-ordering, i.e., A- H>B-H>A-G>... etc. (These arithmetic 
differences are the Improvement Weights. ) 

(3) The preferred or best set was assigned goodness value 100. 

(4) A void was assigned goodness value 0. 

(5) The minimum improvement was made greater than a number 6 > 0, fixed in 
advance, (65 may be chosen arbitrarily, but in practice it is usually chosen between 2 and 6.) 

Kruskal and Aumann of the Princton Analytical Research Group have shown that 
it is possible to obtain the entire space of goodness-value sets satisfying (1) through (5) above, 
A mean point may then be chosen from this space in some suitable manner. 


The Function by Which Over-All Military Worth Is Found 

Improvement Weights having been established, it becomes possible to use these weights 
in arriving at MIP (or ship class-requirement importance) Weights. First, however, one must 
determine the function by which the weights of the various factors will be combined to arrive at 
an over-all value of military worth. 

In order to analyze some of the possible functions, weights were assigned to all factors 
in one specific problem and input data from which transportation cost- matrices could be pre- 
pared were supplied to the David Taylor Model Basin. There a number of different solutions 
to this one problem were obtained by use of the UNIVAC. To obtain these solutions, multi- 
plicative, additive, and combined (multiplicative and additive) formulae were tried. Further- 
more, power and multiplying constants were varied. 

One of the functions tried was 





Worth = MIP Weight + a(Improvement Weight) - (Utilization Weight)”, 


where the constants a and b were varied. 
The multiplicative formula 


Worth = (MIP Weight) . (Improvement Weight)* - (Utilization Weight)? 


also was tried. 


In these solutions the sum of the military worths associated with the allocations is con- 
sidered as the criterion of desirability. : 4 

The solutions prepared by using the UNIVAC were compared with intuitive- method 
solutions prepared by hand in the Bureau of Ships. 

Meanwhile, Weyl] of the Office of Naval Research suggested that the Utilization of 
Features factor might be dropped out, so the problem was brought to the attention of the 
Princeton Analytical Research Group. 

A UNIVAC solution was then run which ignored the Utilization of Features factor. This 
solution corresponded most closely to the intuitive-judgment solutions and seemed to allocate 
equipment in such a way that good utilization of features was consistently made. 





4improvement weights were assigned by methods described previously. MIP weights were 
assigned more or less arbitrarily. Utilization weights were assigned on the basis of the cost of 
features wasted compared with the total cost of the set. 
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Simultaneously, the Princeton Analytical Research Group analyzed all solutions and, on 
the basis of formal theoretical methods, came to the following conclusions: 

(1) Improvement Weights arrived at by the methods previously described should not be 
raised to a power; that is, their exponents should be 1, regardless of the function by which 
weights are combined. 

(2) There are theoretical and practical objections to all functions tried except the 
multiplicative function; most of the faults of the previous solutions were traced to the use of 
the additive formula. (One of the advantages of the multiplicative formula becomes apparent 
when one analyzes the assignment of MIP Weights.) 

(3) Good Utilization of Features is obtained in solutions based on the formula 


Worth = (MIP Weight) - (Improvement Weight) 


and therefore should not be included as a separate factor in the function. 

(4) There are practical and theoretical objections to including the Utilization of 
Features factor in the formula. 

A solution using the formula 


Worth = (MIP Weight) - (Improvement Weight), 


when carefully compared with the hand-prepared intuitive- judgment solution, disclosed no 
objectionable differences that were not directly corrected by changes in the MIP Weights. 
Furthermore, this UNIVAC solution disclosed a number of obvious but previously undiscovered 
mistakes in the hand-prepared solution. At this stage of the development of the problem, it 
became possible to arrive at reasonable methods for assigning MIP Weights. 


The Assignment of MIP Weights 

With Improvement Weights established and a multiplicative formula chosen, it is pos- 
sible to arrive at MIP Weights by means of a few rank-ordering type of decisions. The system 
in use at this time, suggested by Kruskal and Aumann, involves asking an appropriate 
authority questions of the form: 

"If you had one model A set available, would you rather use it to replace a model C in 
the highest priority ship class-requirement or a model F in the lowest priority ship class- 
requirement ?" 

If the authority answers that he would rather replace the model C, then we know that: 





(MIP Weight highest priority) - (Improvement A - C) >(MIP Weight 
lowest priority) - (Improvement A - F) 


MIP Weight highest priority Improvement Weight A - F 
MIP Weight lowest priority ~ Improvement Weight A - C 





where the improvements A - F and A - C correspond to the appropriate ship class- 
requirement. 
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By tabulating all Improvement Weight ratios such as is shown on the right-hand side of 
the above inequality and asking the appropriate questions, the MIP Weight ratios can be obtained 
within narrow limits. By this process the next question might be: 

"Would you rather replace a model B on the highest priority ship class or a model F on 
the lowest priority?" 

If the answer to this question is that replacement of model F is preferable, then we know 
that 


Improvement Weight A- F MIP Weight highest Improvement Weight A - F 
Improvement Weight A- B~ MIP Weight lowest © Improvement Weight A - C 








and the spread of the MIP Weights is established. If the answer is that replacement of B is 
preferable, we might next ask for the preference between replacing model B on the highest 
priority or model G on the lowest priority. 

Likewise, questions can be asked to establish the ratios between the highest MIP pri- 
ority and other priorites within the list. By this method, the entire range of MIP values can be 
obtained. The number of questions can be considerably reduced, however, if it.is possible to 
obtain agreement that a certain range of MIP priorities should all have approximately the 
same weight. (In one problem now being solved, the MIP priority structure was divided into 
six groups. All of the priorities within a group were considered to fall "with but after" the 
highest priority in the group. This reduced to fourteen the number of questions that had to be 
asked. ) 


Vil. EXTENSION INTO THE FIELD OF PROCUREMENT PLANNING 

It can be shown that allocation plans resulting from weights assigned as described above 
are consistent with all the basic rank-ordering decisions. In this sense, the allocation plans 
are valid, and it can be reasoned that insofar as allocation planning is concerned the military 
worth values are valid. It seems reasonable, then, that these values can serve as a basis for 
developing combined "procurement-allocation" plans. This has been done by considering 
(instead of the quantities q; available for allocation from stock or previously placed contracts) 
the quantities G,+ X;, where Xj represents the additional quantity which will be procured. 

The quantities 40, 116, 12, etc., in Figure 1 illustrate the quantities (q's) of sets cur- 
rently on order or in stock. The columns or slots represent an allocation possibility, i.e., a 
requirement-ship class model-installed combination. The quantities y, indicate how many 
such identical allocation possibilities there are in the fleet. Wp}; represents the worth of 
assigning one model D set to slot j. Zn; represents the number of model D sets assigned to 
slot j. 

The procurement-allocation problem can presently be solved in two different ways. 
By one method the amount in dollars available for additional procurement is set as a restric- 
tion, i.e., one restraint becomes 





X1 Pa + Xo Pa +... etc. << amount in dollars available for additional procurement 


where P A is the price in dollars of each model A set, Pp is the price in dollars of each 
model B set, etc. 





tot 
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The solution to this problem pro- 
yides the quantities of the various models 
which should be procured, as well as an 
optimal allocation plan consistent with 
quantities in stock, on order, and to be 
procured. It also is possible to set the 
total military worth desired as a restric- 
tion and obtain a solution which minimizes 
total cost. 

The procurement-allocation prob- 
lem is no longer a type of transportation 
problem, because the quantities available 
are not fixed.> But it retains many of the 
features of the transportation problem. 


Vill. PRESENT STATUS 

One procurement-allocation prob- 
lem has been solved by usingthe analytical 
procedures describedherein. The solution 
was quite satisfactory and differed con- 
siderably from an_intuitive-judgment 
solution previously prepared for the same 
problem. The analytical solution not only 
allocated available stocks in a different 
way butalso used most of the money avail- 
able for procurement to buy different 
models than those bought by the intuitive- 
judgment solution. Furthermore, it pro- 
vided sets for more ships than did the 
latter solution. 
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Figure 1 - Development of procurement-allocation 
plans: (a) allocation, (b) procurement-allocation 


Figure 2 shows the relation of the total optimal military worth sums to the amount in 
dollars required to achieve the worth. Comparison of this curve with others to be prepared for 
different programs which are competing for the limited available procurement funds should aid 
executives in deciding how much money to allocate to each of the programs. 


——— 


SSuzuki has formulated the procurement-allocation problem as a general linear pro- 
gramming problem. 





J. W. SMITH 


MILITARY WORTH OF OPTIMAL 
PROCUREMENT ALLOCATION PLAN 
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Figure 2 - Cost of military worth 
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INTRODUCTION 

Whenever a large quantity of property is under the management of a single organization, 
a well-devised plan for maintaining the quality and keeping accurate account of the property 
can lead to substantial savings. The purpose of this paper is to present several applications of 
mathematics to the solution of these "cycling" problems. The formulations given here are of 
the simplest possible kind. The treatment presented does seem worthwhile, however, for it 
suggests the outlines of solutions to real problems and the feasibility of approaching more 
complicated problems of this type. 

Several problems will be treated. The common element in these problems is that a 
solution consists in choosing a "program," i.e., a sequence of times, at which certain routine 
operations should be performed in order to minimize losses. The solutions presented will all 
be based on the assumption that certain costs functions are known. In some situations we will 
need only to make "geometric," e.g., concavity, assumptions about the costs functions. Inother 
situations we will require the functional form or even the completely specified cost functions. 
In applying the techniques of this paper the acquisition of this information might well be the 
principal difficulty. 

The three main problems considered are concerned with the choice of optimal proce- 
dures in a replacement program, in an inventory control program, and in a preparedness 
program. In the replacement problem the solution involves the minimizing of the long-run 
average-cost-per-unit time of "bank" replacements. The solution to the inventory problem is 
the program which will minimize the costs per year that arise from taking inventories and 
from imperfect information regarding the inventory. In the preparedness problem the solution 
program gives the procedure which will minimize the expected costs required to send equip- 
ment into the field when an emergency arises. 


A REPLACEMENT PROBLEM! 
A decision has been made to replace all the florescent tubes in an office building at one 
time, i.e., a "bank" replacement policy. The problem then is to determine the appropriate 


times to perform these replacements. 


*Manuscript received 2 August 1956 
1A similar problem has beentreated at the Operations Research Groupof the Case Institute 
of Technology. 
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Let us assume that the cost of replacing the tubes (including the costs of the tubes) 
is A. Let F(x) represent the total, "integrated," loss if x units of time are taken between 
replacements. These costs arise from the annoyances and poor working conditions that occur 
when some of the tubes have failed, and from the cost of replacing individual tubes as they fail, 


A program for this problem consists of a sequence of periods Xy,Xq,-- +> Xj, - 


i 
denoted by X, and replacements are made at the times Xj, Xj + Xo --, z Xj» ere 
j= 


loss is measured by the average loss per unit of time and is given formally by 
n n 
C(X) = Limit | 2 (A+ F(x;)) L x | - 
1 1 


n—-o 


We will consider only cases where F(x)/x tends to infinity as x tends to infinity. More 
precisely we exclude the case limit F(x) = minimum At F(x) | 
Xx—»@ * Oo<x<cw 
possibility that any of the x;'S should be chosen equal to infinity. We will also assume that 
F(x) is continuous, non-negative, nondecreasing, and F(0) = 0. 
Under the above assumptions (A + F(x))/x assumes its minimum value. Call this value 
C*, and let S be the set of x's for which the minimum is obtained. Then 


n n 
C(x) 2 ie (A + rap) fe s| 


This restriction precludes the 


n n 
>|2 X; o*|/|E = 
1 / Ll 


unless each x;,...,X, is an element of S. The average costs can be held down to c*, and 
clearly there is no advantage in choosing periods not in S. Thus the problem is solved by 
finding the set S and evaluating C*. If S contains only one element, the solution is unique. 

The fact that all of the periods in the solution need not be of the same length at first 
may be disturbing. One difficulty is that the solution consists simply of a statement of what the 
periods are and not of the order in which they should be used. Actually, however, in practice 
this freedom would be advantageous, since it would give some play in fitting the replacement 
work into other schedules. 


dF(x) 
dx 


If F'(x) = exists, then the critical equation is 


A+ F(x) - x F'(x) =0. 


The derivative (if it exists) of the critical equation is - x F"(x), where F''(x) = 
a2 F(x)/dx?. Thus if F(x) is convex there is a unique solution. The case of F(x) being 
concave is precluded by the condition F(x)/x-s@ as x-so0. 

With this in mind, let us consider several examples. 
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If F(x)= x“, k >1,the critical equation is 
A+ xk -k xk =0 


and the unique solution x* is 


1/k 


Ga) 


Then x* is an increasing function of A and are x* =o. Finally C* = x(;+) 
— 


(1/k-1) 


If F(x) = ef% 4 » 9>0, then the critical equation is 


A+ e%*_1-9xe%- 0 


e%* (9x-1)=A-1, 


which has a unique solution obtainable by numerical analysis. The solution increases with A 
and decreases with @. 

The following is an example where S contains more than one element and thus the 
solution is not unique. Let A= 1 and 


0 0<x<l 
F(x)=< x-1 if 1<x<2 
eX - e241 2<x<o. 


Clearly F(x) is continuous, non-negative, nondecreasing, and F(0) =.0. We see that 


0<x<1 


fa 
x 
1 


if 1<x<2 


2, Xx 
-@ + 
aaa 2<x<o. 





If x¢(1,2), then C(x) > 1, and thus S is the interval (1,2). 

In this section we have assumed that the cost of replacement is precisely A and the 
integrated loss after waiting x units of time from the last replacement is precisely F(x). In 
Practice, the cost of replacement and F(x) will be random variables. We obtain exactly the 
Same solution, however, if we restrict attention to the "expected" loss and now assume that A 
is the expected cost of replacement and that F(x) is the expected cost of waiting x units of 
time, Besides, in the replacement problem the costs are formed by adding together many 
Small costs, and thus by the law of large numbers these costs will be close to their expected 
Values, with high probability. These reasons will be tacitly used for not discussing the random- 


ness of costs in all of the following problems. Figure 1 illustrates the cost function considered 
here, 


419839 0 -57 -4 
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Figure 1 - Illustrative cost functions. 

Each function is infinite for x = 0 because the fixed cost, A, is 
spent instantaneously. The three curves are parabolic in shape, 
but 3 does not have a unique minimum. All the curves go to infinity 
as x does, but the rate of growth of 2 and 3 is more rapid than that 
of 1. The shape of 1 near its minimum is flatter than that of 2, and 
hence precise adherence tothe optimal procedure is not as important 
in 1 as in 2. Curves like 1 could arise when the costs associated 
with a failure are proportional tothe time that the failure (individual 
failures not replaced) is presentand the number of failures per unit 
of time is nearly constant. Cases like 2 couldarise in a contagious 
situation, i.e., where a failure could tend to cause failures in its 
neighbors. Curve 3 was designed to illustrate the situation where 
S contains more than one point. Cases like 3, however, could arise 
where initially the failure rateis very low, then rises a little faster, 
and then finally becomes explosive as in 2. 


An interesting class of problems of the same form as the replacement problem are 
concerned with the optimum times for abandoning a base of operations in favor of a new one. 
In some lumbering situations the procedure is to build a temporary camp in the middle of the 
forest and then to commence cutting the trees that are closest to the camp. In this way a 
circle of ever-increasing radius of felled trees will be made about the base. Eventually, how- 
ever, this operation will be found unprofitable, since excessive amounts of money and time will 
be spent in traveling from the base to the current logging operation and in returning the logs 
to the base. The problem then is just how large a circle should be formed before moving the 
base of operations. 

Another problem2 of this type is met in mining operations. The question there is what 
is a reasonable length for the drifts (tunnels leading away from the shaft) before it becomes 
desirable to dig another shaft. 

In particular, we will discuss a vein deposit of lead ore. It is assumed that the cost of 
digging the shaft and equipping it is $500,000. Further, that for each foot of cross section, 





2Data required to study this problem were supplied by Professor Fred L. Humphrey, ané 
Messrs. Thomas Eisner and Arthur Grenier of the School of Mineral Sciences, Stanford 
University. 
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150 tons of ore can be removed from the vein. Then since the drifts can be dug in both direc- 
tions along the vein and away from the shaft, it is possible to remove 300 tons for each 
increase of a foot in hauling distance. The cost of moving a ton 1000 feet is 1.1 cents. This 
does not include the fixed costs of loading and unloading but does include empty return trips 
and depreciation as well as wages and power. 

With a constant rate of digging and constant number of tons per foot of vein, the time 
length of the operation, the lineal length of the drifts, and the tonnage removed are all propor- 
tional. Thus the minimization of costs per ton removed, per unit time, or per unit length of 
tunnel is just one problem. Therefore, let us consider the problem of minimizing average 
costs divided by length of drift. It should be pointed out that there will be certain fixed costs 
independent of the length of the drift, and all that we are doing here is looking at the variable 
costs. 

The constant A corresponds to the cost of digging a new shaft and equipping it, which 
has been assumed to be $500,000. The function F(x), where x is measured in feet, is obtained 
in the following manner. If the drift is of length x, the cost of extending it to length x +A 
equals the number of tons of ore in a slab A units thick times the cost of removing the ore. 
This quantity is the derivative of F(x) and equals (in units of dollars) 300 A (1.1)(10" 5x), and 
thus F(x) = 1.65(107°)x2, The critical equation is 


5 x 10° + 1.65(107°)x2 - 3,3(107°)x? - 0. 


Solving for x, we find the optimum length of drift to be 1.7 x 104 feet. Clearly, drifts longer 
than 3 miles will not occur in practice (veins would not be good for such distances) and thus the 
problem of excessive costs due to hauling distances is not a limitation in mining operations of 
the type under consideration. 


AN INVENTORY PROBLEM 

A certain warehouse must take a complete inventory of their stock on the first of the 
year, for various fiscal reasons. The management of the warehouse wishes to consider the 
advisability of taking inventory more than once a year. 

Let us assume that the cost of taking an inventory is A. On the other hand, if an 
accurate knowledge of the inventory is not available, certain losses are incurred, e.g., inability 
to fill orders, overstocking, unregulated pilfering, improper rotation of stocks. Let F(x) rep- 
resent the total loss, "integrated," due to inaccurate information regarding the inventory if a 
period of length x is used between inventories. It is reasonable to assume that F(x) is con- 
tinuous, non-negative, nondecreasing, and F(0) = 0. 

A program consists of choosing periods X1>-++X, Such that their sum is 1 (year). 
An inventory is then made at Xj, X,+Xo,---, 1. The total costs of inventories, i.e., due to 
taking them and due to the imperfect information between them, is 


nAt+ re F(x;) ; 
i=] 


The solution to the problem is that program which will minimize this cost. Management has 
available the choice of n and of the periods (x,, varus x): 
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Thus if F(x) = Bx(B >0), the total cost of the inventory program (during one year) vil] 


n 
be nA+ S Bx;. The assumption has been made that an inventory must be taken at the end of 
i= 


n 
the year so that x,=1 (x; > 0) and the loss becomes nA+ B. Thus, when F(x) = Bx, the 
i- 


loss does not depend on the lengths of the periods but only on the number of periods and is an 
increasing function of this number. The solution to this problem is to choose n = 1, i.e., only 
perform the inventory at the end of the year in which case the loss is A+ B. 

In general, if F(1) = B, then we know that the losses may be held to A+ B by taking 
only one inventory at the end of the year. On the other hand, we must lose at least nA if n 
inventories are used. Thus n must be chosen so that 


A+B 
nA<A+B or nsa® =| A |. 


(where [y] is the largest integer <y). The upper bound, n*, for the number of inventories will 
in practice often be excessive. Knowledge of "fk however, provides a systematic approach to 
the solution of the problem. 

In general, the solution to the problem can be found in the following manner. For each 
value of N<n* minimize the expression 


Cy (x) = NA+ F(X) + F(xy9) +. sac F(Xyy) 


N 
subject to the constraints that x,,, >0 for i=1,2,...,N and }) x,,=1. Fora particular 
Ni pas <8 o, *Ni 


value of N this minimization can be carried out with the help of a Lagrange multiplier. One of 

the first things noted is that F'(Xx1) “Peet y (py) at the minimum. This, coupled with the 

constraints, will usually yield the minimum without difficulty. When F'(x) changes sign only 

a finite number of times and is not constantly equal to zero in an interval, then no complications 

will arise. The solution to the problem is then the sequence (Xy4> cs oe Xnn) which yields the 

smallest Cyy(x). In the examples presented here the solutions are all of the form Xn = *y2? 
-+ = Xyy- But examples can be constructed to show that this is not a general result. 

We will now consider two interesting special cases where the solutions are easily 
obtained. If F(x) > Bx, then only the inventory at the end of the year should be taken. This 
condition will be satisfied whenever F(x) is concave. To prove this, we note that for any 
program the cost satisfies 


n n 
nA + Z F(x;)< nA + zL x, B= nA+ BS A+B. 
= i= 


Examples of functions F(x) that satisfy this condition are: 
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F,(x) = 1 - (1-x)¥ k >1 
1 , 
F,(x) = x+ 4 sin (5x) , A>0O 


F(x) = 1 - e 9X @>o. 


If F(x) is convex, the solution consists of N equal periods where N is the integer 


n n 
which minimizes n(a+ (2). The definition of convexity is zu A, F(x;) >F( 2, di x) . 
i= = 


n 
where 0< u< 1 and )) 


: 1 m1 1 
X di = 1. Hence, choosing AG = we have 2 = F(x;) >F( Er = x;). 


i=1" 
n n 1 
Now, using zo x; = 1, we obtain Le F(x,)> n F(-) . Thus the loss when n periods are used 
i= i= 
is minimized when all of the periods are equal. 


In the convex case it can happen that N = 1, e.g., when F(x)= Bx. The minimum can be 


found by evaluating n (a + F( 1 


*) for n= 1, 2,...,n* ina finite number of steps. When F(x) 


is differentiable the critical equation3 is 
x F' (x) = A+ F(x). 


This equation has at most one solution in (0,1), and N must be one of the integers adjacent to 
the reciprocal of the solution. If the equation has no solution, choose N = 1. 
Thus if F(x) = x*, (k >1), the critical equation is 


k k 


kx =A+x , 


. As long as k - 1< A, we have N= 1. 


k-1) 1/k 
and N must be an integer adjacent to (5) 


Other examples for this case are 


F(x) =e -1 9>0 


F9(x) = 1 - cos (5 ax) O>r<1 


F(x) = x“ e% @>0,k>1. 


It is interesting to discuss these functions in terms of the inventory problem. The 
function F(x)= Bx corresponds to the situation where the annoyances due to poor information 
regarding the state of the inventory arise at a constant-rate, no-aging process. The concave 
case occurs when there is a relatively rapid rate of loss at the beginning as compared with that 
at the end. The concave function F,(x) =l-e 6x. for instance, is the probability that an error 





3The same critical equation arose in the Replacement Problem. 
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has been made in recording the stock 
level by the time x when the number of 
stock transactions is large and the prob- 
ability of making an error during a par- 
ticular transaction is small. When there 
are several items in stock the function 


-6 ix 
Zz Bi (1 -e ) could be used to repre- 
sent the expected losses from errors, 
where B; is the cost when the jth stock 
is in error and 6; is the error rate for 
the jth stock. The convex functions cited 
can be rationalized in much the same way 
If F(x) is like curve I, only the annualinven- as the functions used in the replacement 
tory should be used. The relatively rapid in- bl Th 1 hich ith 
crease in the loss at the beginning makes it proolem. © GXAMPlES WHICH AFC NKR 
unprofitable ever to start over again. Curve II convex nor concave correspond to situa- 
is the case that F(x) = Bx, and again only the It tely si d 
annual inventory should be used. If F(x) is like tions where there are a ernately Siow an 
curve III, then it may be desirable to take more rapid increases in the losses because of 
than one inventory. The relatively slow increase : 
in the loss at the beginning might make it advan- poor information about the inventory level. 
tageous to start over again rather than to reach This type of phenomenon could arise 
that portion of the curve where the losses , , 
increase rapidly. where there are periodic fluctuations in 
the demands on the warehouse or if stocks 
are replenished at various fixed times. 
The following problem was suggested by Professor Humphrey. Let us suppose that we 
wish to build a pipeline between two points separated by one unit. The question is how many 
crews should we use on the job. To solve the problem, we must enumerate the costs which 
depend on the number of crews. Let A be the fixed cost of equipping a crew. Assume a crew 
can build a units of the pipeline per unit of time, a crew spends b dollars per unit of time, 
and the "cost" of one dollar for t units of time is ect - 1. Since the pipeline cannot be used 
until it is completed, it is clearly desirable that all crews should begin and end their work at 
the same time. If this is done, the total costs, which are a function of the number of crews, 
are given by 








Figure 2 - An-inventory problem. 


T 
nA+ n | bet - a)at, 
0 


1 
where n = numbers of crews, and T= => is the time required to complete the operation. After 


simplification, this formula reduces to 


n(a+ a ie 


The problem is then solved by finding the n which minimizes this last expression. To a first 
approximation the minimizing n is of the form 


1,/bec 
9A ° 
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The above approximate solution has intuitive appeal, for it says that the number of crews used 
is an increasing function of the interest rate and of the rate of spending money and is a 
decreasing function of the speed of doing work and the cost of establishing crews. 

The analogy between this pipeline problem and the inventory problem is strong. The 
sole important difference is that in the pipeline problem we do not consider the possibility of 
using the crews for periods of different lengths. Interesting modifications of the pipeline prob- 
lem would consider the possibility of using parts of the pipeline as they are finished. 


A PRODUCTION PROBLEM 

This is a generalization of the replacement problem. Its solution is similar to that of 
the former. 

Assume we have two machines to use in a certain part of a continuous production 
process. When the first (second) machine has been in use x(y) units of time we have suffered 
losses F(x)(G(y)) because of the wearing or decalibration of the equipment. Let us assume 
that F(x)(G(y)) is continuous, non-negative, nondecreasing, and F(0) = 0(G(0) = 0). Also 
assume that F(x)/x(G(y)/y) tends to infinity as x(y) tends to infinity. Let A(B) be the cost of 
repairing and of replacing the first (second) machine by the second (first). 

A program now consists of two sequences of numbers, X = (x,, Xo, - - .) and 
Y= (Vv, Yoo: + .). The average loss per unit of time is measured by 


n—>o \i= 


C(X, Y) = Limit (z (A+ Bt F(x;) + ay)/ r (x; + v) ; 
1 i=1 


A solution consists of programs which minimize C(X, Y). As in the replacement problem, we 
see that none of the x; or y; can equal infinity in the solution. The solution is found by 
minimizing 


(A+ B+ F(x) +G(y))/(x+ y), 


and the minimum average loss per unit of time, . is the minimum of the above expression. 

If F(x) is identically equal to G(x), A need not equal B, then the above expression 
becomes (C + F(x) + F(y))/(x+ y), where C = A+ B. Then, from the reasoning employed in the 
replacement problem, it is clear that in the solution we may choose x = y and that the solution 
is found by minimizing (C/2+ F(x))/x. When F(x) and G(y) are differentiable, the critical 
equations are 


(x+y) F'(x) - [A+B+ F(x) + G(y)] = 0 
(x+y) G'(y) - [A+ Bt F(x) + G(y)] = 0 


and the solution to the problem will be found among the solutions to this pair of equations. 

The problems treated in this section are easily generalized to the situation where there 
are several machines instead of just two. The analyses of those problems are exactly the same 
as this one, 
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THE INVENTORY PROBLEM (GENERALIZED) 

We will now present several generalizations of the inventory problem. These gener- 
alizations are much like that made of the replacement problem. We first present the mathe- 
matical problem and then suggest possible applications. 

Let F(x) and G(y) be continuous, non-negative, nondecreasing, and F(0) = G(0) = 0, and 
let A and B be positive constants. Assume that 


n m 
uxt 2k yj=1, 
isl j=1 


where each Xj and yj is positive. Choose the x,'S and y;'s and n in order to minimize 
n m 
nA+mB+ 2) F(x,)+ 2 Gly), 
iz] i=1 


when m= 1, m=n-1, and m=n. 

The case m = 1 could correspond to the inventory problem where thé costs associated 
with the annual inventory are different from those associated with inventories taken during the 
year. The case m =n - 1 suggests a process where we alternate between losses arising from 
F and G and where we want to begin and end in the F process. Thus during a workday it 
might be desirable to alternate working periods and rest periods. Then F(x) could represent 
losses due to fatigue and G(y) could represent the loss due to nonproductivity. It is clear that 
it would be desirable to begin and end the work day with a work period. Finally m =n suggests 
an alternating process where we begin in one phase and end in another. This type of problem 
might arise in therapy, for instance, where we wish to alternate passive and active exercise. 
Assume that the active exercise must follow a period of passive exercise and there is no 
advantage in finishing the treatment with passive exercise.* F(x) might correspond to losses 
in not using the more highly desirable active exercise, and G(x) might correspond to the 
fatigue from the active exercise. 

An analysis of these problems will not be made, since it readily follows from the 
methods already developed. These problems can also be easily generalized to the case where 
there are several processes. 


A PREPAREDNESS PROBLEM 
A certain piece of equipment is kept in storage to be used only in the case of emer- 

gencies, During this period it costs A units of money to inspect the equipment and F(x) units 
of money to repair (like new) the equipment if it has been x units of time since the last inspec- 
tion and repair. The costs A and F(x) apply only if the inspection is made on a routine basis 
and not under emergency conditions. Under emergency conditions the corresponding costs will 
be denoted B and G(x). We will assume that A and B are positive constants and that F(x) and 
G(x) are continuous, nondecreasing, and non-negative. Also F(0) = G(0)= 0. It is also reason- 
able to assume that G(x) > F(x); and in many applications G(x) will be considerably larger than 





*“Historical evidence supporting this assumption is given by Lewis F. Richardson in 
‘The Distribution of Wars in time,’ Journal of the Royal Statistical Society (N. S.) 107, 
(1944) pp. 242-250. Further development of Richardson's ideas is contained in James R. 
Newman's The World of Mathematics, Volume 2, (1956) New York, pp. 1238-1263.” 
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F(x), particularly for large x. Finally, the model is completely specified by saying that the 
probability that an emergencv occurs in the next T units of time is 1 - e 47 where A>0O and 


is known. A program now consists of a sequence of numbers X = (x, Xo, - .). Routine 


n 
inspections are made at the times Ky, Xy + Xq,--- > a Xi where the emergency occurs between 


n n+1 
times ). Xj and }> Xis and emergency inspection is made at time T, the time of the emer- 
1 1 


gency. 
The total costs of inspections and repairs from the beginning of the program until the 


emergency is met is 


n 
Z (A+ F(x,))+ B+ G(x"), 


n 
where x* =T- X;- This is a random variable, since T is a random variable (which also 
i=l 
forces n to be a‘'random variable). The expected value of this quantity is found after simpli- 
fications to be 


x: 
-Ay, itl 
e a. F(x,)+ | AG(z)e" 4? ae . 
i=1 0 


xy 
rx) | AG(z)e" 4? az+ B+ 
0 


i 
where y; = Pa x,. A program which minimizes C(X) is a solution.4 If C, (x) is the con- 
J = 


7. 


ditional expected loss incurred after X4> given that T >X4, we have 


"1 -A -A 
C(X)= j AG(z) e” 4? az+ a(t -e * < (oe F(x,)+ C, (x). 
0 


Because of the nature of the exponential distribution, i.e., P(T >x,+t |T >X,) = P(T >t), it 
is clear that C,(x) does not depend on x,. In fact, C,(X) is given by the same formula as 
C(x), where x,,, now plays the role of x,;. Thus if S are the possible values of x, in the 
Solution, then it is clear that Xo as well as all of the later x;'S must be chosen from S. Thus 
the minimum cost can be obtained by choosing all of the x;'S equal to some particular member 
of S, say x. Then C(X) = C, (x) and 


x 
e 4x (A + F(x))+ | AG(y) a i dy 
0 


C(X) = Bt 
1 a e 4x 








4Minimizing C(X) minimizes the expected cost per emergency and the long run average 
Cost per unit of time. 


419839 O -57-5 
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The members of X can be found by minimizing the above expression for C(X), and a solution 
consists of an infinite sequence of these values. The solution will be unique if S has only one 
member. 


Assume 


" A 
| AG(y) eT Y ay 
0 


exists and that C() is the expected loss if no routine inspections are made, then 


CO 


wo 
AG(y) so dy - e74* (a+ F(x) + | AG(y) e4Y cy) 
0 


x 





C(w) - C(K) = a 


If the numerator of this last quantity is always less than zero, then no routine inspections 
should be performed. ‘ 


If C'(X) = a C(x) exists, the critical equation is 


x 
(1 - ce" 4*) (F'(x) + AG(x)) = A(A+ F(x)) +A j, ae"4Y Gy) ay. 


In order to simplify notation, let us assume that time is measured in the proper units to make 
A= 1. Then, after integration by parts, the critical equation becomes 


x 
(1 - e™*) F'(x)+ G(x) - F(x) { eY G'(y) dy= A. 


Consider the case F(x) = 0, i.e., the equipment can be inspected and repaired for a 
fixed cost if done routinely. Then the critical equation becomes 


x 
G(x) - j, eY G'(y) dy=A. 


This critical equation has at most one root. Thus if G(x) =\x, this last equation becomes? 
AX ~a(1 - e-*) =A. 


The solution can be found by numerical methods. 





>The following equation and equation (5) on page 283 of Dynamic Equipment Policy by George 
Terborgh, published by McGraw-Hill (1949), are of the same form. This suggests a strong 
similarity between the problems discussed by Terborgh and the present preparedness problem. 
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If F(x) = G(x) the critical equation is 


x 
(1 - e”*) F'(x) - { eY F'(y) dy= A. 
0 


The derivative of the left-hand side is 
(1 - e°*) F"(x), 


so that if F(x) is either concave or convex there is at most one solution. Thus if F(x) = 1 - 
eX where A > 0, the critical eauation becomes 


na te . 9) . far (i - e"@+)) x) =a 


2 
-Ax | _.-x a) |. .. 
° le (2 +A} = A+iTy» 


which can be solved by numerical methods. 

In general, the derivative of the left-hand side of the critical equation is proportional to 
(G'(x) + F"(x) - F'(x)). If this function is of constant sign, then the critical equation has at most 
one root and S has one member. This condition will be met if G'(x) > F'(x) and F"(x) is 
positive, i.e., F(x) is convex. 


FURTHER PROBLEMS 

1. One of the most unrealistic features of the models treated inthis paper is that there 
isno "memory."" The functions F(x) and G(x) have been treated as if they do not change as 
the total elapsed time of the program increases or as history of the process changes. Thus it 


n 
might be desirable to replace F(x), and similarly for G(x), by F(x, ~ x,) or F(x, Xj, Rea x)» 


where x is the time since last (x,) inventory, replacement, etc. 


2. The stochastic nature of the problems might be further emphasized. Thus, in the 
inventory problem, instead of making the decision of invéntorizing or not at each instant of 
time, we might consider the possibility of making a decision of whether or not to take a statis- 
tical survey at each instant of time and then, as a result of the survey, decide whether or not an 
inventory should be made. 

So far we have considered only expected losses, but in many situations the variability 
of the losses, which could be characterized by the variance, also might be of importance. 

3. In both the inventory problem and the replacement problem the results should be 
integrated with the existing inventory theory which covers other aspects of the problems. 

4. In this paper little effort has been made to present a theory on how the functions 
F(x) and G(x) and the other parameters should be specified. All that has been done is to give 
plausible forms for these functions. Undoubtedly, for some of the problems, theoretical con- 
Siderations could delimit the possible choices for these functions. In order to apply this work, 
empirical studies, however, would be required. 


‘ 
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The author discusses a simplified mathematical model of the 
allowance list, and draws some general conclusions. 











1, INTRODUCTION 

One of the more difficult problems of naval logistics is the preparation of adequate 
allowance lists for naval vessels.2 On the one hand, one must have reasonable assurance that 
the requirements for consumables and technical spares will be met under most circumstances. 
On the other hand, one is confronted with the severe limitations of available space aboard ship, 
with budgetary considerations, as well as with a host of other less important constraints. 

This problem is normally solved by the application of judgment based on past experi- 
ence. The question arises whether or not it is possible to construct a simple mathemat- 
ical theory of the allowance list problem which will make better use of the available usage data 
related to some definite program elements (activities upon which consumption of commodities 
depend). 

The problem we shall consider can be stated as follows: How does one properly stock 
a vessel with the necessary commodities, subject to the limitations of available space? This 
entails a decision on the number of commodities to be carried, the quantity of each commodity 
to be stocked, and the relative weight (i.e., "military worth", etc.) to to given each commodity, 
all subject to the overriding considerations of space. In this paper we shall attempt to treat a 
somewhat idealized version of this problem by a method which might be extended to more com- 
plicated situations. 








*Manuscript received 5 January 1956 

1Research performed under contract with the Office of Naval Research. The author wishes 
to express her appreciation to Dr. C. B. Tompkins and Dr. M. A. Woodbury for stimulating and 
informative discussions. She is also indebted to Dr.I. Heller and Dr. W. H. Marlow for reading 
the manuscript and for their valuable criticisms. 

2By an allowance list we mean a listing of distinct commodities (including all classes of 
naval material, regardless of cognizance) which should be aboard ship for the maintenance of 
the ship for a specified time period. This listing contains also the quantity of each of these 
commodities to be stocked. Each commodity is uniquely identified by a specific Standard Navy 
Stock Number. 
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We must now define an objective which we shall attempt to achieve by our construction 
of the allowance list. If we knew exactly what the demand for each commodity would be, we 
would have no problem. We would simply stock the required amounts aboard ship, if space 
permitted, or partially aboard ship and partially aboard supply ships or shore facilities and 
deliver to the ship as its supplies diminished. 

The fact of the matter is that we cannot predict the demand with certainty. We can, 
however, estimate the probability distribution of the demand for each commodity for a given 
type of ship, under certain known operating conditions, from a statistical analysis of usage data, 
We can use this information to construct our allowance list in such a way as to minimize the 
probability of depletion of any one item in a given period. This might be a reasonable criterion 
for commodities of critical importance. Alternatively, we can choose another criterion to 
guide us in the construction of our allowance lists. In this paper we shall attempt to construct 
our allowance list in such a manner as to maximize the average number of demands fulfilled. 
This is an arbitrary criterion at best. It is a reasonable choice, however, for those items 
whose usage varies fairly widely and which are not to be replenished by a supply ship at fre- 
quent intervals. 











2. STATEMENT OF METHOD 

Let us outline our method. Consider a given commodity, say a. Let the amount stocked 
be yy and the amount demanded be x,. Let ¢, (xq) be the probability that an amount x, will 
be demanded during the period under consideration. For the sake of simplicity, we shall con- 
sider the case where x, is a continuous variable and ¢q(Xq) has the analytic properties 
required for the existence of derivatives and integrals of the function ¢,(xq). Then the average 
amount demanded will be: 


© 
E(x,) = | Xqy by (Xq) dx,- 
0 


If the demand x, is less than y,, this demand will be met. But if the demand x, exceeds the 
amount stocked (y,,), then only an amount y, can be met; a demand of (x, - y,) is then left 
unfulfilled. On the average, the amount supplied (when y, is stocked) is: 


Ya rr) 
(2.1) Uy (Yq) = j Xqy Pq (Kg) dXq t Ya by (Xq) dXq >» 
0 Yo 


and the demand left unfulfilled is, on the average, 


(2.2) Vay (Yq) = [ (Xq - Yq) bq (Kg) dXq = E(Xq) - Ug (Yq). 
Ya 


We assume the demand for commodity @ to be independent of the demand for commod- 
ity 8. Thus, the average amount supplied of commodity @ is dependent only on yg, not on 
anything pertaining to commodity 8 . In practice, this assumption may not always be correct, 
but it will greatly simplify our analysis. 





AN OPTIMUM ALLOWANCE LIST MODEL 


Our problem now is to choose the set Yy» Yor --- Y, 80 as to satisfy the condition 
r c 
Ca. Va = 
“= 


(i.e., the total cube of items stocked must equal the available space), and simultaneously to 
maximize 


U(y,, eee > Y,) = yr Uy (Yq) » 
a=1 


which is the total quantity supplied on the average. Here, c, is the cube occupied by one unit 
of commodity a; and C is the total cube available forthe storage of all the n commodities. In 
general, we may wish to attach a different weight, w,, to each commodity, depending on its 
importance (i.e., "military worth," etc.) and its unit of issue. We do this by defining 

Wiy,, ii yy) as follows: 


n 
Wy, oe 9 yy) a Zz Wa Ug l¥ gq) 
a=1 


and maximizing W, subject to the condition 


We shall refer to Wy, peels yy) as the "total utility function;" and we shall call wgUg (Yq) 
the "@ -th utility function." 

Once we have maximized Wy, Dales y,) subject to the space constraint, we know the 
individual Yq'S. Given these y,'s, we can compute the average number of the demands for 
each @ that we can fulfill, namely Ug (Yq)'S, aS well as the average number of the unfulfilled 
demands, the vq (yq). From the ug(yq) we obtain the optimum value of Wy, er Yn)» which 
we shall denote by Wo» and call it the "maximum utility function."" It should be clear that Wo 
is a function of C only, since for a given C the Yq'S, and therefore the u, (y,)'s, are fixed by 
the maximization procedure. The significance of Wo will be elucidated in the discussion (see 
Section 5). 


3. RESULTING FORMULAE 
The details of the calculation will be relegated to Appendix I. The resulting formulae 
for the Yq'S are: 


Voy Ca 
j o (Xq) dx,=1-’—; a 
0 Va 


w 
~ Cy Yo= C, where min (ee 
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These are (n + 1) equations in the (n+ 1) variables Yy> Yor +++ 9 Vy and A (where A isa 
Lagrange multiplier). They determine the amounts, y a 5 which are to be stocked. Using 
these values of y,'s, we can compute u, (y,) and v, (y,), for each a, from equations (2.1) 
and (2.2). 

For some simple distributions, like the exponential, these equations can be solved 
explicitly for the y,'s and’. For the more interesting distributions, however, the explicit 
solution of the equations is not feasible, and an iterative procedure is indicated. We shall 
illustrate the meaning of these equations and the method of their solution by working out sev- 
eral examples. 

In the first example we shall assume all the $, (x,)'s to be normal. This is a fairly 
important case, since the normal distribution is the limiting continuous form of a wide class 
of distributions which occur in practice (e.g., the Poisson distribution). 

The second example will illustrate the properties of this model when all the ¢q (xq)'s 
are logarithmic normal distributions. The reason for this choice is that most of the tests of 
modern statistics are based on the assumption that some function g(x) of the variable x is 
normally distributed. In this case we take g(x) = In x. This distribution is easy to handle 
analytically, and the numerical results provide an enlightening illustration of the principles 
involved. 

Finally, we shall consider a case where some of the ¢, (x,)'s are normal while others 
are logarithmic normal. This will show how the allowance list is to be constructed when com- 
modities are distributed according to several different distributions. 


4. NUMERICAL EXAMPLES 
Example 1 - The Normal Distribution 


Let ¢q(Xq), for reasons stated above, be the normal probability density functions. 
Then equations (3.1) become: 


Ce 1 (Yq-Mg)/oy 


2 
1-A See et /2 at 
Qa 


-0O 


n 
LX caYa=C;, 
a=l 


where m, and a2 are the mean and the variance, respectively. 


In order to solve equations (4.1), we choose a value of in the permissible range, 
compute the y,'s, and from their values determine the corresponding C. In practice, how- 
ever, we are given C, not A. We therefore repeat this calculation for several values of \ until 
we find a A which corresponds to the given C. We are then ready to determine u,(y,)'s, 
namely the average amounts supplied. These can be written as: 


2 2 
fa Co Sq -(Ya-MQ)"/20 
netted me -rgzoralt- Ga sbae 
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We shall be able to gain further insight into the meaning of these equations for yq and 
Ug(¥q) for the normal case if we assign definite numerical values to the parameters of these 
distributions and compute the resulting y, and uy, (Yq). Such a numerical model is presented 
in Table 1.1. It should be noted that the mean is taken to be rather large, since we are con- 
sidering the normal distribution as the asymptotic form of some other distribution. It is only 
in this case (i.e., when m,/o q@ > 4) that equations (4.1) and (4.2) are sufficiently accurate for 
our purposes, 


TABLE 1.1 
2 E (xq) Cq = 2400 





E (xy) = mg | VVar (xy) =O, 


i) 
R 





100 3 
100 10 
100 3 
100 10 
100 3 
100 10 
100 3 
100 10 


onroawr won 
NNN NY SY  S 
ane Ke oo 























These numbers were picked so as to represent quantities with a wide range of com- 
binations of cube, "weight," and variance. In order to acquire a feeling for the order of magni- 
tude of a given C, one should compare it with 2 E (xy) cg. This sum would bethe C if one 

a 


were to take yg equal to the mean usage of commodity @, for all @. In Table 1.2 we list the 
results, ? 


TABLE 1.2 





d= 0.01; A= 0.02; d= 0.05; d= 0.15; 
C = 2698.3 C = 2647.5 C = 2565.4 C = 2409.7 





Ya | Ya Yq) | Ya | “ala)] Ya | “a Vo Yo | & OS 





99.99 99.98 99.94 
99.97 99.93 99.79 
99.94 99.86 99.55 
99.79 99.53 98.51 
99.99 99.99 99.97 
99.98 99.97 99.91 
99.97 99.94 99.81 
99.91 99.79 99.38 


on ow Fr wwe 









































419839 O -57 -6 
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Because the E(x,) = 100 for all a inthis example, u, (y,) also is numerically equal 
to the percentage of demands met. This also holds for examples 2 and 3, with the exception of 
the first six cases of example 2. In that case, E(x,) = 10, and therefore 10u, (y,) is numer- 
ically equal to the percentage of demands met. We shall postpone the detailed interpretation 
of these numerical results to Section 5. 


TABLE 1.3 





AW 
AC 





2409.7 | 1180.7 


2481.5 | 1189.8 


2565.4 | 1195.9 


2647.5 | 1198.7 





2698.3 | 1199.4 














Example 2 - The Logarithmic Normal Distribution 
Let us consider the case where 1n Xq is normally distributed.3 Then 


1 -(1n xq - m_)2/20q" 


2 
=————- e here and o@ are E(1nx_) and 
dy (Xq) tno, x, ,w My (1n x,) 


a 
Var (ln x,,), respectively; and equations (3.1) become: 


1 (iny,- m,)/o 
a a’ a 72/2 ce 
or e cee oe pee 


- © Wa 


’ 


Le =C. 
= a Ya 
Also, u, (y,), the total quantity of commodity @ supplied on the average, is given by: 


(In Ya - Mq - Oy /8y t2/2 1 cO 
u = E(x,) — e dt + ——7 
a Va) = E(%q) fe va | 
_ 


% 


2 
et /2 at, 


m 


My +02 /2 


E(x,) =e 





3See, for example, A. Hald: Statistical Theory with Engineering Application, p. 160. 
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In Table 2.1 we shall illustrate this model by means of a numerical example. It should 
be noted that this numerical model differs from the one given in Table 1.1 in including cases 


with a much smaller mean. The results appear in Tables 2.2 and 2.3. 


TABLE 2.1 


= E (xq) Cq = 1980 
a 





= 
fe) 


E (x,) NVar (x,) m 0 


Qa a 


R 
R 





10 3 2.25949 | .29359 
10 10 1.95600 | .83256 
10 30 1.15129 | 1.51742 
10 3 2.25949 | .29359 
10 10 1.95600 | .83256 
10 30 1.15129 | 1.51742 

30 4.56208 | .29359 
100 4.25860 | .83256 
300 3.45388 | 1.51742 

30 4.56208 | .29359 
100 4.25860 | .83256 
300 3.45388 | 1.51742 


—_ 

KF OoODMDDA DOL WD & 
ee ee ee ee 
annrre KB OOO eS ee 


_ 
i) 





























TABLE 2.2 





d= 0.01; d= 0.02; r= 0.05; d = 0.10; rd = 0.15; 
C = 6261.18 C = 4522.98 C = 2706.96 C = 1712.63 C = 1318.47 





Tn iM DE Ta i Me Bell Se i te Melt te ie eT em | 





9.98 9.91 9.89 9.78 9.66 
9.81 9.67 9.31 8.79 8.32 
8.84 8.47 7.43 6.28 5.43 
9.89 9.78 9.40 8.63 7.52 
9.31 8.79 7.47 5.56 3.68 
7.43 6.28 4.20 2.23 0.99 
99.79 99.10 98.93 97.80 96.60 
98.14 96.71 93.08 87.89 83.20 
88.44 84.69 74,26 62.80 54.37 
98.93 97.80 94.03 86.34 75.25 
93.08 87.89 74.71 55.61 36.84 
74.26 62.80 41.96 22.27 9.93 


oon nowt won 


—_ — 
no KF © 
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TABLE 2.3 





Wo 





1318.5 


1712.6 


2707.0 


4523.0 





6261.2 





391.79 


453.98 


524.67 


581.89 


607.90 











Example 3 - The Case of Several Distributions 
In this model some of the $, (x,) are normally distributed, while others have the 
logarithmic normal distribution. In Table 3.1, the first 6 cases have the logarithmic normal 
distribution and the last4thenormal. y,, and u, (y,) are computed as in Example 1 for 
= 1,..., 6 and as in Example 2 for a =7,..., 10. The results are tabulated in Tables 3.2 
and 3.3. 
TABLE 3.1 


= E(x,) c,, = 3000 
a 





= 
R 


- 


E (x,) 


N Var (x J 


My 


o 
a 








ooaosnsouwrtrwown 


— 
oO 





le 





aqqarredgdgaa&#aeaw = = 





100 
100 
100 
100 
100 
100 
100 
100 
100 
100 





3 
10 
30 

3 
10 
30 

3 
10 

3 
10 





4.60472 
4.60020 
4.56208 
4.60472 
4.60020 
4.56208 

100 

100 

100 

100 





0.03000 
0.09994 
0.29359 
0.03000 
0.09994 
0.29359 
3 

10 
3 

10 
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TABLE 3.2 





— 


A= 0.01; 
C = 3647.30 


A= 0.02; 
C = 3496.25 


r = 0.05; 
C = 3268.53 


d= 0.10; 
C = 3049.40 


r = 0.15; 
C = 2860.67 





Va 


Ug (Yq) 


Ya 


Uy (Yq) 


Ya 


Ug (Yq) 


Va 


Ug (Yq) 


Va 


Ug (Yq) 





oon out wn 


_ 
o 


99.99 
99.96 
99.79 
99.93 
99.75 
98.93 
99.99 
99.97 
99.94 
99.79 


99.98 
99.91 
99.58 
99.85 
99.44 
98.04 
99.98 
99.93 
99.86 
99.53 


105.01 
117.28 
155.25 
102.00 
106.39 
116.76 
104.94 
116.45 
102.02 
106.75 


99.93 
99.75 
98.93 
99.54 
98.32 
94.02 
99.94 
99.79 
99.55 
98.51 


99.85 
99.44 
98.04 
98.78 
95.82 
86.34 
99.86 
99.53 
98.80 
96.01 


99.76 
99.13 
96.59 
99.41 
91.70 
75.58 
99.77 
99.22 
97.53 






































91.76 





TABLE 3.3 





Wo 





2860.67 


3049.40 


3268.53 


3496.25 














3647.30 





5. DISCUSSION 

Some fairly general rules can be gleaned from a scrutiny of these tables. It becomes 
clear for the normal case that commodities with a large variance are stocked to a considerable 
degree in a large vessel, but only very sparingly when the available space is rather limited. 
In fact, when the space becomes very hard to get, one stocks considerably less than the amount 
equal to the mean usage for those commodities which have wide fluctuations of demand. In 
Appendix I we show that a modified version of this rule has a rather wide range of applicability. 
We show also that for many distributions, when the mean is multiplied by k and the variance by 
Ke then yq and Ug (yq) also are multiplied by k. This also is noticeable in the numerical 
model (see Table 2.2). It also can be seen that commodities with large cube must be stocked 
Sparingly in a small vessel. These results are true for a wide class of distributions. The 
Specific numerical values of the y,'s and the Uy (¥,)'S are, of course, more sensitive to the 
particular distribution and may have to be computed for each individual case. 
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The very fact that the general rules which come out of our model are eminently plau- 
sible serves as an a posteriori justification of the basic assumptions we have made in the 
Introduction. 

We come now to the discussion of the behavior of the maximum utility function, Wo: 
We note that as C (total available cube) is increased, Wo at first rises sharply, but then its 
rate of increase diminishes. This becomes especially clear when one examines the behavior 
of the ratio AW/)/ACc. The behavior of Wo can be used as a guide in deciding on the amount of 
space C to be set aside on a given ship for storage of commodities. 

From the above discussion it appears that the present model leads to reasonable results 
for commodities of a noncritical nature. It might, therefore, be used as a basis for some pre- 
liminary allowance list computations. 

It is clear, of course, that our criterion of maximizing the average number of fulfilled 
demands may not be the only useful criterion. For example, one might use the criterion that 
the sum of the mean square deviations of the demands from the amounts stocked shall be a 
minimum, subject to the space constraint. In other words, we could minimize 


2 ™ 2 
E E (Xq - Yq) | af | (Xq - Ya) Pq (Kq) dXq 
0 


= 2 {Var (X,) + [E(xq)-Yal?| 
a 


subject to 


This would yield the relations 


n 
C- E 
PP CB (x. ) 


Yq = E(x,)* cy - or @ o4, ...68. 


2 
£1 8 





c 


Because this model minimizes the mean square deviation of the demand from the amount 
stocked, it could be expected to lead to overemphasis on meeting the demands of commodities 
with very large variance at the expense of those commodities whose demand is reasonably 
uniform. Our model, on the other hand, does not minimize this mean squared deviation, but 
concentrates instead on meeting the greatest possible number of demands on the average. 
Thus, it appears that the criterion that we have chosen, while not unique, does lead to desirable 
consequences, 

The model that we have considered in this paper is a purely probabilistic one. We have 
assumed that the probability density functions were known. In practice, the probability distri- 
bution is not known, and only some of its properties can be obtained from the analysis of the 
available data. One can use this statistical information in order to determine the parameters 
in an appropriately chosen distribution function. This, in fact, is the way in which one must use 
the present model for detailed computations. On the other hand, one could conceivably 
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reformulate the problem so as to take more direct advantage of the available information, pos- 
sibly bypassing some of the intermediate steps inherent in the present formulation. 


APPENDIX I 
The General Theory 
Let us define the following set of functions 


Wy Xy When xy4<y 


a “a a 


fy (Xa » Yq) = a=1,...,m, 
Wa Yq When Xy 2 Yq 


where x, is the amount of commodity a demanded, and y, is the amount of commodity a 
stocked, and where wy arethe relative weights assigned the commodities according to their 
importance. This function, fy (xq, Yq), represents the "gain" that we achieve if we meet a 
demand x, when we stock an amount y,. Of course, when x, is larger than y,, we meet the 
demand as best we can, namely, by supplying what we have in stock yy. 

Our objective is to maximize our total gain. We approach this by maximizing the total 
utility function, W (defined as the expected value of the total gain), subject to the space con- 
Straint. In other words, we maximize 


n 
W (y;, +--+» ¥,) = =|, {, @ v0) subject to z 9, *% . 


It should be noted that we have assumed the gain to be additive. Now, 


E|= f_ Xqs ¥a)| = ff Z fy (Xqs Yq) © (Xj, --- » X,) AXy,--- dX, 


where ® (x), nee x.) is the joint probability density function of Xo +++ 5X): 
We have made the assumption that the demand for commodity @ is independent of com- 


modity 8 , so that their joint probability density function can be written as: 
® (x), see y x.) =o, (x1)o5 (x5) - O45 (x,) . 


where $(Xq) is the probability density function of x,. The range of each xq Clearly is from 
0 tom. Whenever no limits of integration are indicated, they are to be taken as from 0 to o. 
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Since [ea (xq) dxq =1, 


we get: 
E[E by Oy + Ve)] = ty Oxy » ¥y)y Oy) dry +... + fy Oy» Yq) Oy OQ) dX, - 


Now let us define 
Wa Yq (Yq) = E fa Ka Va) = [fa Kas Ya) a (Ka) AXq 
and therefore maximize 
Wy, 9 ows y= 2 Wy Uy (Ya) subject to Bla ¥q=C- 
Using the method of Lagrange multipliers, we form the function 
D= 2 Wq Yq Ya) ->[E ea Ya -C]- 
a a 
We get n equations of the form: 
aD 
) 


d 
11 e—_ -—_—— - e 
(I. 1) 0 Ye a¥e [Wy Yq (Yq) - X&g 


The set of equations (I.1), together with the equation of constraint, determines the n +1 vari- 
ables y, and’. Equations (1.1) can be written as: 


d 
Wa dyq ° (Ya) 


d Yo i>) 
* aye Xa be (Xy) dx, +Y, { oy a) oa] 
0 Yo 


@ 
{ oo (x,) dx, , 
Vo 


Remembering the normalization of the by (x,), we have the final form of our equations: 


Ya t. 
(1.2) | %o (x,) dx, z1-A— 


0 Wa 


=C. 


and zc 
a a Yo 


These equations determine a unique set of Yq'S, provided that, for all a,¢, (x,) > 0 for all 
values of xq (except possibly on a set of measure zero). 
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That this stationary point is really a maximum of Wy), re yy can be easily shown 
by the methods outlined in Courant-Hilbert, 'Methoden der Mathematischen Physik," 1931 
edition, p. 200. 

We are now ready to compute the mean quantity of commodity a supplied when an 
amount yg is stocked. This is obtained as follows: If an amount of x, is demanded which is 
less than y,, then that amount is supplied. If, on the other hand, the amount demanded, X,» is 
larger than the amount stocked, only an amount Yq can be supplied. Thus the average amount 
supplied is: 


a rs) 
Uy (¥ g) f Xy bog (Xq) dXq + Va bq (Xq) AXq . 
0 Yo 


The mean quantity of commodity @ not supplied when an amount y, is stocked is: 


Vay (Yq) = [ (Ky - Yo) by (Xq) dXq - 
Voy 
Clearly, 
Uy (Yq) +Vq (Yq) = E(Xq) 
where 
E (x,) = |raba (x,) dx, 5 


In order to solve equations (I.2), we shall assume a value of 4; compute y, by using 
a table of the appropriate cumulative distribution function; and then determine the corresponding 
C from the subsidiary condition. Since, 


Vo 
ox] ox, dx, <1, 


Ww w 
it follows that = >rA >0 for all a, and therefore min( —*)> 4 >0. This defines the range 
of permissible values of ). 

We shall now turn to the discussion of several useful properties of our model. We shall 
Show that if the mean is multiplied by k and the variance by k“, then Yq and u a also are 
multiplied by k. This is usually true for distributions which have more than one parameter. 
For example, it holds for the normal and the logarithmic normal but does not hold for the 
Poisson. The proof proceeds as follows (for simplicity we delete the subscript a): 

Given a distribution function ¢(x), such that 


(1) [oo ax- 1 
(2) E(x) = [xo (x) dx= M 


(3) Var (x) = Es (x) dx - [E(x]? = 2 a 
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then the distribution function 


W(x) = 54(Z) ,  k>0 


has the properties: 


(1) fv (x) dx = feo) dx = [ore az- 1 
(2) E(x)= [xvi ax- {xzo(Z) dx 


= k [2 ¢(2)dz= kM 


1 /x 
(3) Var (x) = [ev (x) dx - [E(x)]* = fx? xo(;) ax - [E(x)]? 
= k? [Pew dz - [kM]? =k? >?. 
If the distribution function (x) has the same analytical form as ¢ (x), differing from 


it only in the values of its parameters, then we can use it to find the scaling laws for y and 
u(y). Thus, if y is defined by the equation 


y c 
I @ (x) dx = 1-Az, 


y' y' 7k 
roe fe aaa xo (i) ox-f ¢(z)dz, 


and therefore 


This defines the new amounts to be stocked. Also, if originally 


y co 
u(y) =f xo(x)ax ry [ o (x) dx. 
0 y 


' 


' y - 
u(y) = | xw(x)dx+y' { w (x) dx, 
0 y 


f HG 


0 


J k (ve) 
u'y’ =k fi soGdés+y' [ o(z) dz= ku(y). 
0 y /k 
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It also is easy to show another fairly general property of the model. Ifa function g(x) 
of the variable x is normally distributed, with mean m and variance o”, then 


g(x) - m 
SS 


is N(0, 1). Then, 


2 2 
— -(g(x)-m)"/20" dg (x) 
1-7; * —~ 


We shall consider the case where g(x) is a monotonically increasing function of x such that 
-0 <g(x)<@, while O0<x<o. Then 


~ 2 
f, (x) dx == — aie es /2 dé=1- n= 
-@ 


c 1 
Clearly, whenever 1 - weg we have g(y) =m, which can be solved for y. If 1-A + >3 ; 
f — 
then giy-m) = T >0, and therefore g(y) = m+ to. Clearly 7 is the same for two commodities 


having equal = , even though their means and variances may differ. But if T, = To, My = Mo, 


and 0; >», then clearly g, (y,) >85 (Yo); and therefore also y; >Yo- Similarly, if 


c U1 gly)-m 
1- A, <3, Gp = - T<0; and if T,= Tg, My = Mg, and 0; >Op9, then 81 (¥y) <Bq (Yo) 


and ¥y<Yo- This shows that for small C (large positive X) one stocks less of the commod- 
ities with large o. On the other hand, when C is large, one stocks more of the commodities 
with large o than of those with small o. 








THE BASIC THEOREMS OF REAL LINEAR EQUATIONS, INEQUALITIES, 
LINEAR PROGRAMMING, AND GAME THEORY 


David Gale* 
Brown University 


1, PREFACE 

The purpose of these notes is to present short and direct proofs of five important 
theorems of real linear algebra, three existence theorems concerning solutions of linear 
equations and inequalities, and the basic theorems of linear programming and game theory. 
Needless to say, none of these theorems is new, and in the case of the existence theorems it 
would probably be impossible to determine when they first became known. In the last decade, 
however, they have assumed a new importance because of their application to a wide class of 
problems in fields connected with mathematical economics, engineering, logistics, etc. Since 
the same theorems occur over and over again in these applications, it seems desirable to 
record them in an organized manner for reference purposes. 

It should be emphasized that we have made no effort to give "best possible" proofs 
from a pedagogical point of view, because a number of expositions of this sort already exist. 
Thus, our treatment will not contain any numerical examples, plausibility arguments, dia- 
grams, bibliographical references, or convex cones. This presentation is intended rather for 
the reader who is willing to take his algebra straight. 

On the other hand, an effort has been made to make the presentation as direct as pos- 
sible. There are no "prerequisites" except familiarity with the algebraic properties of the 
real numbers. There are no preliminary lemmas, and no concepts are introduced but those 
directly involved in the theorems. 

As to the choice of theorems, the inclusion of the results on games and programming 
requires no justification. Concerning the existence theorems, a word should be said. 

Let us consider the problem of determining whether or not a set of simultaneous 
equations has a solution. In case the system does have a solution, the proof of existence can 
be carried out "constructively" by the simple means of exhibiting a set of numbers and verify- 
ing that they satisfy the equations. On the other hand, if the system has no solution, there is 
no constructive way directly from the definition to prove nonexistence, for it is clearly impos- 
sible to examine in turn each of the infinite number of candidates for a solution and show that 
it fails to satisfy the equations. Now, the existence theorems in the next section are intended 
to remedy this situation by giving constructive criteria for ascertaining when a system of 
equations or inequalities has no solution. The results actually take the form of nonexistence 
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theorems, which say that a given system has no solution if, and only if, a certain "dual" system 
has a solution. 

A similar situation obtains in the case of extremum problems. In order to show from 
the definition alone that a given function attains a certain value , as its maximum, one would 
have to examine all other values of the function and verify that they were no greater than p. 
This impossible task is circumvented in linear programming by the duality thvorem which 
asserts that is the maximum of a certain function if, and only if, it is the minimum of the 
"dual" function. Now the problem can be solved constructively by simply exhibiting the argu- 
ments for which the values of the two functions are equal. (Note, on the other hand, that the 
problem of showing that a given value of a function is not a maximum can be solved construc- 
tively by simply exhibiting a larger value of the function.) 

In the next section we list the familiar vector and matrix notations to be used in the 
statements and proofs of the theorems. While this may seem unnecessary, there is still so 
much diversity of notations, even in such a standard subject as linear equations, that once 
again for the sake of completeness and convenient reference we include our own glossary. The 
reader is assumed to be familiar with these matters and is encouraged to skip the next section 
and get on to the theorems, turning back if necessary to the section on definitions if some step 
in a proof is not clear. 

The theorems of the last three sections are listed in logical order, each depending on 
its predecessor, except for Theorem 2, which is not needed for the results that follow but is 
included because of its importance in many other applications. We have given each theorem 
a title in order to stress the logical sequence of ideas. 

Finally, we remark that the only property of the real numbers used in our proofs is the 
fact that they form an ordered field. Thus all our theorems are also valid, with the real num- 
bers replaced, for instance, by the rational numbers. As an application of this fact it follows, 
for example, that a game with a rational payoff matrix has a rational value. 


2. DEFINITIONS AND NOTATION 

An n-vector x is an ordered n-tuple of real numbers (E4, itn Ey) The number bi 
is called the ith coordinate of x. We shall use the symbol 0 to denote the n-vector all of 
whose coordinates are zero. The symbol e, the ith unit vector, is the n-vector (21, nee é.) 
such that & i? 1, £, = 0 for j4# i. Note that a 1-vector is simply a real number. We shall 
consistently denote vectors by Latin letters and real numbers by Greek letters. 

The set of all n-vectors is called real n-space, denoted by R,- We define three 
operations on n-vectors. 


(1) Addition. If x = (2;,---5&,) and y= (nj,... sn)s then 








x+y= (4 +1y>--- »E,+n,)- 


(2) Scalar Multiplication. If x = (24, onwlg ll n) and A is a real number, then 





AK= (AE},...,2E,). 





(3) Scalar Product. If x = (E4, preree Ey) and y = (ny, --- »n n) are vectors, then 


n 
(x, y) = Zz Ein. 
i=1 
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We shall take for granted all the properties of these operations which follow imme- 
diately from the definitions and will use them freely in the proofs of the next sections without 
giving explicit justification. Examples of such properties are: 


(A +u)x=Ax+ux, 
A(x+y)=AX+AY, 


(x, Ay)= (Ax, y) = A(x, y) y 


(x, ¥y + Yo)= (x, yy) + (%, Yo), 


(x, x) 20, and (x, x) = 0 if, and only if, x = 0. 


We also define a relation between n-vectors. 

(4) Partial Order. If x= (é),..., &) and y= (nj,... »1,), then x<y and y>x if 
ESny at oc so 

One immediately verifies among others the following properties: 

If x>y and y>z, then x>z, 


If x; >¥, and X5>Yo, then x; +X9>¥,+Yo, 


If x>y, then Ax>aAy for ,;>0 
Ax<Ay for <0, 
If x>y and a>0O, then (a, x) >(a, y). 


We shall fly in the face of convention and call the vector x positive if x >0. The usual 
term for such vectors is "non-negative.'' We have chosen the unusual terminology to avoid the 
frequent use of this awkward double negation in what follows. Note that in our terminology the 
vector as well as the number zero is positive (and also negative). 

An mxn-matrix is an "array" of numbers is i=1,...,m,j=1,...,m. The ith 
row vector a; of A is the n-vector (a5) weeny @),)- The jth column vector a’ of A is the 
m-vector @, gore @ ). The number a ij is called the ijth coordinate of A. 

The transpose A’ of the matrix A is an nxm-matrix whose ijth coordinate is a 
where or =a 








’ 
ij’ 
ji 

(5) Matrix Multiplication. 

If A is an mxn-matrix with column vectors al saetae a” and x = (4, eee e,) is an 
n-vector, then 





Immediate properties: 
If ai is the ith row vector of A, then the ith coordinate of Ax is (a;, x). 
A(x+ y)=Ax+Ay, 
A(Ax)=AAx, 


and especially important is the following: 
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If x is an m-vector, y an n-vector, then 
(Ax, y) = (x, A’ y). 


3. EXISTENCE THEOREMS 
In this section, A will be an mxn-matrix whose ith row and jth column are, respec- 
tively, a; and al: x is an m-vector; and y an n-vector. 


THEOREM 1: (Positive Solutions of Linear Equations) 
The equation 


Ay=a 
has no positive solution if, and only if, there exists x such that 
A’ x<0 and (a, x)>0. 
PROOF: If y is a positive solution of (1) and A' x <0, then 
(a, x) = (Ay, x) = (y, A’ x)<(y, 0)= 0, 


so (2) has no solution. 
The converse is proved by induction on n. For n= 1, (1) and (2) take the form 


(1); 1 =a 


(2), (al, x)<0 and (a, x)>0. 


1 1 


Assuming (1), has no positive solutions, if a° = 0 then x =a satisfies (2),. If a 70, 
let X, = (a?, a) a- (a?, a) al and observe that (al, Xy) = 0. If (a, x,) 7 0, then one easily 
verifies that x = (a, x1) X, Satisfies (2),. Finally, if (a, X,) = 0, then (x,, X,) = (al, a!)(a,x,)- 
(al, a) (al, X1) = 0 so x, = 0, and hence a = ((a?, a)/(a}, a) al Since we have assumed (1), 
has no positive solution it follows that (al, a) <0, hence again x= a satisfies (2),. 


Now assume the theorem true for a matrix of n - 1 columns. For convenience we 


write a = a? 


n-1 
. By the hypothesis of the theorem the equation )° 1 j al = a? has no positive 
j=1 


solution, hence by the induction hypothesis there exists X, Such that (ad, X,) <0 for j<n-1 
and (a°, x1) >0. If also (a", X,) <0, then x= x, satisfies (2). If (a", x1) >0, define 


al. (a", X,) al - (ad, x1) a" for j=0,....n-1. 
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n-1 
a - 
j=1 


(1') 


has a positive solution, then 


-1 . -1 : 
a? = ‘> nN al ~ 1/(a", x,) z a (a), x1) aa (2°, x)| a" 
j=1 j=1 


is a positive solution of (1), contrary to hypothesis. Therefore, applying the induction 
hypothesis to (1'), there exists Xo such that 


(2') @!, x.) <0 and (@°, x») >0 OP Pe Bice OSE. 
Let x = (a, X1) Xo -@, Xo) Xy and verify 
(al, x) = @, x.) for j = 0, coe y D- 1, 


(a", x)= 0. 
Then (2'), (3), and (4) imply (2), completing the proof. 


THEOREM 2: (Solutions of Linear Inequalities) 
The inequality 


Ay<a 
has no solution if, and only if, there exists a positive x such that 
A'x=0 and (a, x)<0. 
PROOF: If y is a solution of (5) and x is positive and satisfies A'x = 0, then 
(a, x) >(Ay, x) = (y, A'x)= (y, 0) =0, 
so (6) has no solution. 


Conversely, suppose (6) has no positive solution. That is, letting a= (a peers an)? the 
equations 


(5') 


have no positive solution, hence by Theorem 1 there exists a vector Vy and number n such that 
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(a,, yy) + n@,<0,i=1,...,m and - 1 >0. 
Let y= - 1/n yj, and verify from (6') that y satisfies (5). 


THEOREM 3: (Positive Solutions of Linear Inequalities) 
The inequality 


Ay<a 
has no positive.solution if, and only if, there exists a positive x such that 
A'x >0 and (a, x) >0 
PROOF: If y is a solution of (7) and x is positive and satisfies A'x = 0, then 
(a, x) >(Ay, x) = (y, A'x) >(y, 0)= 0, 


so (8) has no solution. 
Conversely, if (7) has no solution, then the equation 


m 
(7") , 14 ¥ av, e-a 


j=-1 ? i=1 


has no solution y, >0, A; >0 where e, are the unit vectors in R,,. By Theorem 1 there 
exist x, such that 


(8') (aj, X,) <0, (e;, X,) <0 and (a, X}) >0. 
Let x= - X,, and verify from (8') that x satisfies (8). 

4. DUALITY THEOREM OF LINEAR PROGRAMMING 
Let A be an mxn-matrix, a an m-vector, b an n-vector. 
Let Y be the set of all positive solutions of 


(9) 


and let 


Let X be the set of all positive solutions of 
(10) A'x >b, 


and let v = min (a, x). 
x 
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THEOREM: A necessary and sufficient condition for the existence of 
or v is that both X and Y be non-empty, in which case yp = v. 


PROOF: We first show the necessity of the condition. If X is empty, then clearly v 
does not exist, so suppose x, «X and Y is empty. Then (9) has no positive solution, so by 
Theorem 3 there exists Xo >0 such that A’ Xo >0 and (a, Xo) <0. Hence, for any positive 
\, X= X, +AXy satisfies (10) and (a, x) can be made arbitrarily small by choosing ) suffi- 
ciently large, hence v does not exist. 

A symmetric argument gives the corresponding result for yp. 

Assuming now that X and Y are non-empty we prove the existence and equality of pu 
and v. First observe that 


(11) (b, y)<(a, x) for all xeX,yeY, 


for, 
(b, y)<(A'x, y) = (x, Ay)< (x, a), 


so it remains only to show that for some xe X and ye Y, (a, x) <(b, y), which is equivalent 
to the statement that the inequalities 


(12) Ay<a, - A'’x<-b, (-b, y) + (a, x)<0 


have a positive solution. If they do not, then by Theorem 3, there exists a positive m-vector u 
and n-vector v and a number n >0 such that, 


(13) A'u - nb>0, - Av+na>0O and (a, u) - (b, v)< 0. 

First, suppose n = 0. Then, from (13), A'u >0 and - Av >0; but since (9) and (10) 
have solutions, it follows from Theorem 3 that (a, u) >0 and (-b, v) >0, contradicting the 
third inequality of (13). 

On the other hand, if » >0 and 1 = 1/nu, V= 1/nv, then from (13) 


A'U>b, AV<a and (a, 0) <(b, ¥), 


but then U eX and VeY and this contradicts (11). Thus (13) has no positive solutions, and the 
proof is complete. 


5. FUNDAMENTAL THEOREM OF GAME THEORY m 
Let X be the set of all positive m-vectors x = (6, vow bn) such that )> Fr =1, 
i=1 


n 
Let Y be the set of all positive n-vectors y =(ny, veeg Ny) such that a ny = B. 
j= 


Let A be an mxn-matrix. 
THEOREM: There exists X, € X, Yok Y and a number © such that 


(x,» Ay) >w>(A'x, Y,) for all xeX, yeY. 
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PROOF: We first consider the case A >0. Let e be the m-vector, f the n-vector, 
all of whose coordinates are 1. Now the set Y; of all positive solutions of Ay <e is clearly 
non-empty (let y = 0) and bounded, for since A> 0, no coordinate of y can exceed ‘ae (a, ip 


Therefore, the numbers (f, y), y « Y; are bounded above, and hence by the duality ashlee 
there exists y,©Y, and x, positive such that A' x, >f and (e, X1)= p = (f, y,) >0. Let 


w= 1/p, X= WX, Yo= WY. Then x, eX, y, e¥, Ay < we, A’x, >wf, and hence if x «X, 
ye Y. Then 


x, Ay) # (A' Xo» y) >(wf, y) =W: (x, we) 


> (x, Ay,) = (A' x, Yo) . 


If A has some negative coordinates, choose v< min a. and let A be the matrix with 


ij 


i, 


coordinates aj = @j;-. Then A >0, so there exists x, ©X, y, e¥, and @ such that (x, Ay)> 


@ >(A'x, y,) for all x eX, yeY. But Ay = Ay - vf and A’ x = Ax - ve, so that letting 
w = @ + v one again obtains (14). 





AN EXTENSION OF JOHNSON'S RESULTS ON JOB LOT SCHEDULING! 
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Johnson's method for efficiently solving certain two-stage production 
scheduling problems is generalized to include cases in which some jobs 
require only a single stage and in which the two-stage jobs may require 
the machines in both of the possible orders. 











An efficient computational procedure for solving the following problem was obtained by 
S. M. Johnson [2]: 

Jobs 1, 2,...,N are to be processed by Machines 1 and 2. Each job n requires a 

single continuous interval of Machine 1 time, of length an >0, followed by a single 

continuous interval of Machine 2 time, of length b, >0. A job can be processed by, at 
most, one machine at a time, and a machine can work on, at most, one job at a time. 

It is desired to assign the required intervals in such a way as to minimize the over-all 

time required to complete all the work. 

Johnson's paper also generalized his procedure to include a very special three-machine case, 
and the present writer has provided a slight additional extension [1]. 

The purpose of the present note is to point out a simple but substantial generalization 
of all these results. We shall present this generalization in terms of a two-machine problem 
parallel to that described above. The reader should have no difficulty in extending it to three- 
machine cases analogous to those discussed in [1] and [2]. 

The generalization consists of removing the requirement that each job go first to 
Machine 1 and then to Machine 2, and replacing it by the demand that each job be of one of the 
following four types: 

A. Jobs which require processing by Machine 1 only. 

B. Jobs which require processing by Machine 2 only. 

C. Jobs which require processing first by Machine 1 and then by Machine 2. 

D. Jobs which require processing first by Machine 2 and then by Machine 1. 

The sets of jobs is still indexed over n= 1, 2,...,.N; and Job n requires an interval of 
length ay > 0 on Machine 1 and an interval of length b, >0 on Machine 2 (a, >0 if Job n 
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actually requires Machine 1, ai: 0 otherwise; and similarly for db, and Machine 2). It is 
again desired to assign the required intervals in such a way as to minimize the over-all time 
needed, 

Types A, B, C, and D include all possible one- and two-operation jobs, such that the 
order of performance of the operations on each two-operation job is pre-designated. 

As indicated in Johnson's paper, our task can be reduced to that of choosing for each 
machine the order in which it is to perform the operations required of it. With this in mind, 
we shall refer to a pair of orderings, one for the collection of operations to be performed by 
each machine, as a"'schedule."" A solution to our problem will be called a "best schedule." 
There is no difficulty in going from a schedule to the best fully detailed assignment of intervals 
consistent with the orderings of the schedule (see [2]). 

It is a simple matter to show that successive interchanges of the positions of jobs will 
lead from any given schedule to one at least as good, relative to our optimizer, and having the 
following properties: 

1, On Machine 1, all jobs of Type C precede all jobs of Type A, which in turn precede 

all jobs of Type D. 

. On Machine 2, all jobs of Type D precede all jobs of Type B, which in turn precede 
all jobs of Type C. ; 

3. The jobs of Type C are in the same order on both machines. 

4. The jobs of Type D are in the same order on both machines. 

We omit the almost trivial proofs. 

Next, it is easy to see that for schedules having the above properties, the orderings of 
the jobs of Type A and B are irrelevant; that is, they have no effect on our optimizer. Thus, 
we are left with the need to choose orderings of the jobs of Type C and the jobs of Type D. A 
proof essentially identical with Johnson's, and hence omitted here, shows that: 

5. There is a best schedule among those with Properties 1-4, and with the added prop- 

erty that the jobs of Type C are ordered by the following rule: Find the smallest a, 
or b, for any job of Type C. If it is an an put the job in first place among the jobs 
of Type C. If it is a b_, put the job in last place among the jobs of Type C. Repeat 
this process until all iSbs of Type C are ordered, at each stage considering only 
those jobs not yet placed and only those positions in the ordering which are not yet 
occupied. 

. There is a best schedule among those with Properties 1-5, and with the added prop- 
erty that the jobs of Type D are ordered according to the rule of (5), but with "first 
place" and "last place" interchanzed. 

Now we summarize the above rema) <s in a computing routine: First, order the jobs of 
Type C and D by the rules of (5) and (6), re .pectively. Then order all operations required of 
Machine 1 by placing the jobs of Type C first, and in the order previously obtained for them; 
the jobs of Type A next, and in any order; and the jobs of Type D last, and in the order pre- 
viously obtained for them. Finally, order all the operations required of Machine 2 by placing 
the jobs of Type D first, and in the order previously obtained for them; the jobs of Type B next, 
and in any order; and the jobs of Type C last, and in the order previously obtained for them. 

It is worth mentioning that when several of the Types A, B, C, D actually appear in the 
problem, there will often be a very large class of optimal schedules. Often a (necessarily 
optimal) schedule according to which each machine works continuously until all the work 
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required of it is done will be obtained by using any order in which Machine 1 does jobs of 
Type C first, then jobs of Type A, then jobs of Type D; and Machine 2 does jobs of Type D 
first, then jobs of Type B, then jobs of Type C. 
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This paper presents tables which give the operating characteristic 
and average sample time functions of sequential probability ratio tests 
concerning the mean of a Poisson process and analogous discrete time 
results for the exponential distribution. 











1, INTRODUCTION AND SUMMARY 

In [1] the study of sequential tests of hypotheses about stochastic processes with con- 
tinuous time parameter was initiated. Among other results, formulas were derived for the OC 
(operating characteristic) and AST (average sample time) functions of sequential probability 
ratio tests concerning the mean of a Poisson process. This process is a familiar one and 


arises in many practical situations. It is generally assumed to obtain when one studies meson 
incidence on a Geiger counter or such phenomena as the initiation of telephone calls and the 
expiration of light bulbs in life testing, or when one counts bacteria in milk on a slide ("time" 
is replaced by one coordinate of the slide). The purpose of the present paper is to give tables 
which will give the OC and AST functions of a wide variety of such tests. 

The OC and AST functions computed from these tables will also prove useful as approx- 
imations to the OC and ASN (average sample number) functions for the now classical Wald 
(discrete time) sequential probability ratio test (see [2]) for testing hypotheses concerning the 
mean of a sequence of independent observations with common Poisson distribution. (See [1] for 
a discussion of the relationship between the classical tests of sequential analysis and the con- 
tinuous time tests of [1].) The tables also yield exact results for discrete time problems 
involving the exponential distribution. 


2. THE TESTS 

Let X(t)(t 20) be a Poisson process with stationary independent increments and mean 
occurrence time 1/A (i.e., EX(1) =A is the expected number of occurrences per unit time2), 
So that the probability that X(T) = k is e"AT (A T)E Ac! for k=0,1, 2,... . It is desired to 
test the hypothesis Hy:A =Ay against Ho:A = )* Ay where Ay >0, *>1. Let Ay, Xo be 





lUnder contract with the Office of Naval Research. 
2We mention the following typographical errors on p.259 of [1]: the mean occurrence time 
was statedas A insteadof 1/A, anda minus sign was printedas aplus sign in equation (4.2) of [1]. 


*Manuscript received 5 March 1956 
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positive numbers with a, +@_ <1. The optimal test procedure (in the sense described in [1}) 

which satisfies the condition that the probability of a wrong decision be a i (or even Sa;) when 

H; is true is described by two real numbers, J >0 and K< 0, as follows: Stop the first instant 
that 


r 
(2.1) ze) -—2 
r 


is< K or DJ, accepting H, or Ho, respectively. (We may take X(t) to be continuous on the 
right, and thus such a first instant to exist, with probability one.) Here 


» . log a* 


(2.2) r “3 me 


K is given by 


=. 
(2.3) K = (to 23 )hona, 


and J is a number satisfying 


1-a5 
(2.4) I< (log ‘- ) logA*< J +1. 
= , = 


The value of J may be obtained exactly as the solution of either of the equations 
(2.5) p(r* ; J, K) = l-a, 


or 
(2.6) p(1/a* ; J, K)= a9 


(substituting herein for A*, K, and @ Or %» as prescribed above). Here the function p (which 
is tabulated below) is given by putting (for A > 0) 


= _ 1 } 1 
(2.7) os (log x)/Q - 1), r + 


r(i) <1, 
(2.8) f(r(A); J, K) =p; J, K), 


and, for r >0, 


= [-(-i)re"T] /i ! 
3 <2 


> [-(-K-i)re"™] / i! 
i<J-K 





f(r; J, K)= e8?. 
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If the tables which follow are to be used to obtain J, a good check is to use one of the equations 
(2.5) or (2.6) to obtain J and then to make sure the other is satisfied. For certain parameter 
yalues one equation may be superior to the other for interpolating in the tables, as will be 
clear from inspection, and this will dictate which equation to use to obtain J. 

It will often be the case that the user of the test does not care about the exact values 
of a, anda » but only about their order of magnitude. In this event one can quickly scan the 
tables of p and find (tabled) values of J and K which give values of a 1 and M» (read off the 
tables at once as indicated by (2.5) and(2.6)) which are reasonable for his use, avoiding inter- 
polation entirely. 

The OC function L(\) = probability of accepting Ho when EX(1) =) is obtained for 
the above test from the same function p (and its tables) through the formula 


(2.10) 1 - L(A)= f(ar*/d4; J, K) = p(q(ar*/,); J, K), 


where q is the inverse function of r; i.e., q(A r*/4) (written as q for short in the equation 
below) is defined by 


Ar* logq 
ee 
1 q-1 


We may also write (2.10) as 


Ay log 


2.11 ,J, K)=1-L 
(2.11) P(q ) es 


This will usually be the most rapid and convenient form for sketching an OC curve from the 


tables. 
The AST function, A(A) = expected duration of test when EX(1)= ), is given by 


A(a)= > h(aQar*/ay); J, K) 
1 
: s(r r*/r, : J, K), 
where h, defined by (2.12), is tabulated below, and s is defined for positive r by 
s(r; J, K)=1(J + 1) - 1(J - K+ 1) f(r; J, K) 


+ er(J-1) } e Ti i (-r)) {_-Kr 
i<J-K-1 jsi 


. £(r; J, K) (J-K-i-1) - m(J-i-1, j)}/i!, 


where I(t) is the greatest integer<t and m(y, j) = y) if y >0 and = 0 if y<0, j being an 
integer 2 0. Corresponding to (2.11), the most convenient form for sketching the AST curve is 
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(2.14) 





(5) F@y 


r* (q-1) h(q; J, K) . 


Ay log q 


3. THE TABLES 

The functions p and h defined in Section 2 are tabulated for values of \, J, and K, 
which, it is hoped, will be useful in many applications. The values of J and -K chosen for 
tabulation were 2, 2.5, 3, 3.5, 4, 5, 6, 8, and 10. Not all possible pairs of (J, K) values in this 
range are tabulated, but rather those pairs where the ratio of the larger to the smaller of the 
numbers J and -K is not too large (usually <2). The tabulated range will thus include (for 
values of \* included in the table) such pairs of values of a 1 and M» as .01 and .05 (see (2.3) 
and (2.4)), but perhaps not .002 and .05 because the last pair would be less likely to be spec- 
ified in an application. 

For each pair (J, K), the values of A for which p and h were tabulated were chosen 
according to the following criteria: (1) only values of A between .05 and 20 were included 
(many practical cases will be covered by this range); (2) for fixed (J, K), the a; being per- 
mitted to range from .05 to .001 (subject to the remark of the previous paragraph), whatever 
values of \* correspond to such a yp % for a given J, K, the tabulated values of \ should 
permit locating three points on the OC and AST curves between the values Ayr Ad: These 
criteria have, of course, resulted in the tabulation of values of A which vary greatly with J 
and K. For example, for J = 2, K= -2, the A values range from .05 to 20, with .4, 1.0, and 
2.0 being the consecutive values which include 1; whereas for J = 10, K = -10, all eleven of the 
tabulated A values lie between .5 and 2.0. The successive values of \, J, and K were chosen 
with the object of making the jumps between successive p and h values regular and not too 
large. But the application of standard interpolation methods to these tables will not be too 
difficult, even though successive intervals between tabled values of A, J, and K will often not 
be equal. As was remarked in Section 2, many users will be able to select an appropriate pair 
(J, K) (giving reasonable values of @ vp a») from a quick scan of the tables, without using inter- 
polation. 

Thirteen figures were used in the computations? (an amount necessary in some cases 
to obtain four or five significant figures for p and h), and the tables are accurate to within one 
unit in the last digit printed. As checks on the computation, each value of p was checked 
against the bounds (easily obtained as a consequence of (2.3), (2.4), and (2.5)) 


74 J+1_, 


A 
ate: J. “esr 
er eo) - S741. 


at #3, 


J J+1 
T-K LP; J, K) <pRG 


Similarly, using Wald's equation (see [2] and equation (2.1) of [1]), it is easy to derive the 
following inequalities on h, which were used to check every computation of that function: 





3The computations were performed onthe Cornell Computing Center's CPC by Mr. Richard 
Lesser and Mrs. Dorothy Hartman, to whom the authors express their thanks. 
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J - (3 - K) pQ3 3, K) fs) 0 y, (J+ 1)-(3+1-K)pQ;J,K) . y 
1 - 1/r (a) {3 h(a; J, K) {3| 1-1) if va {5} 1, 





3? + (K? - 3%) p(1; J, K)<h(1; J, K) 


< (+ 1)? + [K? - (0+ 1)?] p(t; J, K). 


It is interesting to note that the left side of each of the inequalities of (3.1) and (3.2) is 
the classical approximation of Wald which "neglects excess over the boundaries" and which is 
of common use in applications of sequential analysis (to discrete parameter problems). These 
approximations can be rather poor (in terms of relative error) for certain values of the param- 
eters in the present continuous parameter problem. As one would expect, the approximation 
of the left sice of (3.1) is worst when A and J are small, the percentage of excess over the 
boundary being large with high probability in such cases. For large A, however, the relative 
error in the approximation for 1 - p (which is pertinent; see (2.5)) may be even greater. For 
example, if the test J = 2= -K is used when A* = 10, the actual ay is .0046, whereas the 
approximation gives .010. When Ho is true in this last example, the approximation under- 
estimates A(X) by 20 percent. We need not recount the economic losses which can result 
from under- or over-estimation of Ay, Xo and A(A). 

On the other hand, for J = 10 = -K, the relative error of the approximation for a, is 
greatest when \* = 2, a, = .0008, and the approximation gives .0010, while the relative error 
in the approximation for A(X) never exceeds 4 percent. This indicates that the approximations 
will be close enough in most practical applications requiring values of J and -K beyond those 
tabulated. If values of J and -K smaller than 2 are required, (2.9) and (2.13) can be computed 
directly in fairly short time. 

The tables canalso be used to obtain exact results for discrete time sequential problems 
involving the exponential distribution; see Eq. (3). 
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14,000 
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J= 3.5 


K = -2.5 








7.0000 
5.0000 
3.5000 
2.5000 
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1.7000 
1.0000 
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.050000 
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J=5 K = -10 
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NEWS AND MEMORANDA 


Readers are invited to submit to the Managing Editor items of general 
interest in the field of logistics. 


CAPT. C. Stein, Jr., S.C., USN, CDR. M. I. Rosenberg, USN, and Miss J. S. Campbell 
have been appointed Associate Editors. 


The Logistics Research Division, formally attached to the U. S. Naval Supply Research 
and Development Facility at Bayonne, New Jersey, has been incorporated into the Advanced 
Supply System, Research and Development Division (Code SS/W), BuSandA. 


Dr. Merrill M. Flood has been appointed Associate Director of the Engineering Research 
Institute of the University of Michigan and is in charge of the Willow Run Laboratories. A large 
percentage of the Laboratory's activities are connected with Project MICHIGAN, a battle area 
surveillance systems project sponsored by the three armed services and supervised by the 
Signal Corps. An Operations Research Group, established early in 1955 under the supervision 
of D. H. Wilson, is engaged in several governmental and industrial projects in systems engi- 
neering and operations research. Studies of traffic flow and control are included among the 
many other research programs in progress at the Laboratories. 





An International Conference on Operational Research, which is being organized by the 
Operations Research Society of America, the Institute of Management Sciences, U.S.A., and the 
Operational Research Society of Great Britain, will be held at Oxford, England, from September 
2 to September 6, 1957. The object of the conference will be to unify and extend the science of 
Operational Research. To this end, papers emphasizing the underlying unity of Operational 
Research, its methodology, and applications will be presented in four sessions. 


Further information regarding this conference will be circulated later this year. In the 
meantime, proposed papers are invited. Anyone interested should send a summary of his paper 
(not more than 200 words) to Dr. Thornton Page, 7100 Connecticut Avenue, Chevy Chase, Mary- 
land, 
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Opportunities will be made available for people coming from abroad to visit British 
Operational Research establishments engaged in studying specific problems of using specific 
techniques. Arrangements are being made to publish the conference proceedings. 





THE ARMY LOGISTICS RESEARCH PROJECT 


The Army Logistics Research Project has been established at the George Washington 
University, under the direction of Mr. Robert C. Faylor, in response to guidance provided by 
the Army Deputy Chief of Staff for Logistics. The Project was initiated with the aim of pro- 
viding a scientific approach to the design of a modern army supply system, hence the name 
Project MASS. Under Project MASS it was possible to develop, from guidance on future strat- 
egy and tactics, principles upon which a logistics system should be built. From a careful study 
of these principles, the broad characteristics which describe a modern army Supply system were 
set forth as shown below. 


The system should be: 
(a) Responsive to the needs of the combat commander 
(b) Flexible to accommodate mobility 


(c) Dispersed to minimize losses from enemy action 
(d) Efficient to ease the drain on national resources 

(e) Effective to insure accomplishment of mission 

(f) Compatible with concept of joint operations 

(g) Expandable for support of small or large operations 


After the development of this broad expression of characteristics, the next step was to 
divide the total system into logical component parts which can be scientifically attacked and 
solved. Any military supply system can logically be separated into the following three broad 
functions: 


Supply operations, which include supply policy, warehousing, depot operation, 
supply procedures, and organization. 

Materiel movement, which includes all activities pertaining to the movement 
of supplies from the supplier to the consumer. 

Information processing, which includes communications as well as data 
processing and accounting functions. In effect, it is the flow of paper work. 


Project MASS is being developed in three phases. The first phase has been completed. 
It included the development of the supply system concept and the design of a system to accom- 
plish the supply function. This design is in very broad terms and is consistent with the charac- 
teristics of a modern army in combat. Phase two is a test phase in which the design will be 
tested by supplying one commodity to the Seventh U. S. Army in Germany. The commodity of 
Spare parts was chosen, since it presented the greatest challenge in terms of information 
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processing with the least burden in materiel movement. Also, the consumption of repair parts 
during peace time is of the same order of magnitude as consumption in combat. The third 
phase of Project MASS is to theoretically design the complete modern army wae syste 

the basis of available guidance, analysis of future warfare, and information Beta “4 
two. re end result can es modern army supply system capable of providing all aie “ 
~ -iaeteaaa army in the field, operating under the strategic and tactical doctrine of 








RECENT PUBLICATIONS 


GLOBAL LOGISTICS AND STRATEGY, 1940-1943. By Richard M. Leighton and 
Robert W. Croakley, Office of the Chief of Military History, Department of the Army; Super- 
intendent of Documents, Government Printing Office, Washington, D. C. 


This new work in the series entitled THE UNITED STATES ARMY IN WORLD WAR II 
describes the task of effecting the orderly assembly, movement, and delivery of great masses 
of men and materiel throughout the world to meet not only U.S. requirements but also those of 
other Allied nations. 


The authors attempt to analyze in detail the manner in which limitations of supply and 
transportation shaped the strategy of the Anglo-American high command. The work also 
reveals the intricate relationship of logistics and strategy in the top-level direction of the war, 
showing how logistical limitations developed, the degree to which they affected strategic plan- 
ning, and how a variety of such problems were solved. 


The policy, method, and organization which provided the Army's logistical support, 
including the war-time activities of the Army Service Forces, are among the subjects dis- 
cussed. In addition, there is a comprehensive account of American military lend-lease oper- 
ations, including the air program to Russia. 


Included in the volume is the agreement made by the United States to provide enough 
shipping to Britain to continue the war, when it appeared that German submarines would 
threaten her merchant fleet with extinction and the British population with starvation. It was 
this agreement, one of the little-known episodes of World War II, which precipitated a crisis in 
war plans and almost disrupted the Allied high command. 


To many observers, the net achievement of a year and a half of American participation 
in the war seemed discouragingly meager—with fewer forces overseas than during a compa- 
rable period in World War I, dispersed among half a dozen theaters all over the world, and 
with no effective concentration anywhere against enemies who had the advantage of interior 
lines, During this period the Army was mobilizing, deploying, and gaining positions for the 
great offensives that were to come in 1944 and 1945. 
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The United States was confronted during World War II with a new logistical problem, 
one that was interconnected and global, the writers point out. Every local logistical problem 
was part of a larger whole; none could be settled without considering the impact they would 
have on other local problems, often reaching around the world. 


THE QUARTERMASTER CORPS: OPERATIONS IN THE WAR AGAINST JAPAN. By 
Alvin P. Stauffer, Office of the Chief of Military History, Department of the Army, for sale by 
the Superintendent of Documents, Government Printing Office, Washington 25, D. C., 358 pp. 


To win wars, fighting men need not only fire power for combat action and protective 
devices for saving lives, they need also the everyday necessities—food, clothing, petroleum 
products, and items of general utility—without which martial life cannot be sustained. The 
Army's Quartermaster Corps has as its main responsibility the provision of these four classes 
of supply. How it supplied the men who fought under MacArthur and the other commanders in 
the Pacific is the subject of this volume in the expanding series U.S. ARMY IN WORLD WAR IL 


The author describes the vigorous efforts made by quartermasters in the Philippines 
in 1941 to accumulate the food, clothing, and gasoline that would be needed in case of attack 
and shows how these efforts were largely frustrated by lack of time, by the necessity of help- 
ing to supply the 100,000-man Philippine Army, and by the initial American strategy that 
required dispersion of stocks among several depots, which could not be fully evacuated in the 
confusion of the hurried withdrawal to Bataan. He depicts the desperate plight of the defenders 
of that beleaguered peninsula as their food supplies dwindled. He tells of the ingenious efforts 
to stave off starvation by fishing, by harvesting the local rice crop, by slaughtering carabao, 
and most of all, by heroic but finally tragic attempts to break through the strangling Japanese 
blockade and bring in food from the southern Philippines, Australia, and the East Indies. Dr. 
Stauffer explains how, in the end, starvation forced the surrender of Bataan. 


He gives considerable attention to the war's early months in food-importing Hawaii, 
where the Army, fearing a Japanese invasion, gave the Quartermaster, Hawaiian Department, 
an extraordinary role as controller of civilian food. Stress is laid on the importance of Hawaii, 
New Zealand, and particularly Australia as the only large land masses that were developed 
sufficiently to serve as major bases. Since shipping shortages and remoteness from the United 
States put a premium on local procurement of military necessities, Australia and New Zealand, 
as large food producers, proved especially useful. Of the estimated 3,617,000 measurement 
tons of military supplies that Australia had furnished to U.S. armed forces by mid-1945, nearly 
half consisted of food. Without this food American troops below the equator could not have been 
sustained properly. 


The volume covers the supply lines that spread from the depots in the United States to 
widely scattered island bases, most of them 5,000 to 7,500 miles from San Francisco. It 
describes how from Hawaii to primitive New Guinea the difficulties imposed by lack of the 
common storage and distribution facilities available in economically developed lands were 
mostly overcome and supplies brought to troops dispersed over tiny atolls and islands over- 
run with jungles. The occasional interruptions in the flow of military necessities are not 
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ignored but are examined in considerable detail. The author explains the role of Quarter- 
master troop units in keeping front-line organizations supplied, and he discusses the utility 

of individual items under combat conditions and shows why it was sometimes necessary to 
modify them to meet the extraordinary demands of tropical and island warfare. Some parts of 
the volume shed much light on the conditions under which the GI passed his days. Chapter XI 
describes his food and clothing and Chapter IX the provision of baths and other necessities of 


living. 


Dr. Stauffer has delved deeply into the voluminous records relative to Quartermaster 
activities, whether they were conducted at higher headquarters, in bases, or by Quartermaster 
troop units in support of combat operations. His narrative naturally reflects the viewpoint of 
the troops and the commanders in the field. It is one of the very few comprehensive accounts 
of supply operations in a major war. Moreover, it is an interesting account that helps the 
reader understand why logistics is the indispensable companion of strategy and tactics and 
particularly why the satisfactory supply of Quartermaster items is essential to an army's 
well-being. No comparable work on Quartermaster operations in an entire theater of oper- 
ations has been published. 


The volume is one of a four-volume group that records the experiences of the Quarter- 
master Corps in World War II. The first two volumes describe the activities of that service 
in the United States, and another, still in preparation, will tell of its activities in Europe in the 
war against Germany. 


SEA WAR - THE STORY OF THE U.S. MERCHANT MARINE IN WORLD WAR II by 
Felix Riesenberg, Jr., Rinehart & Company, Inc., New York. 


SEA WAR is a vigorous account of the terrible price in blood and ships paid by the 
American Merchant Marine in its World War II effort to maintain the logistics pipeline in the 
face of relentless enemy action. Chronicler Riesenberg, the son of the renowned Merchant 
Mariner, uses the medium of interviews with survivors to place in perspective the motivations 
and attitudes of the American merchant seaman in the light of his sociological and economic 
background. 


TABLES OF THE CUMULATIVE BINOMIAL PROBABILITY DISTRIBUTION, The Staff 
of the Computation Laboratory, Harvard University Press, 1955. 


n 
Let e(n, r, p) = cr (1-p)"-*p™. The function E(n,r,p)= >. e(n,x,p) is the proba- 
x=r 


bility that at least r successes will result in n independent trials if p is the probability of 
Success in one trial. The tables give five place values (rounded) of E (n, r, p) for r = 0(1)n, 
n = 1(1)50(2)100(10)200(20)500(50)1000 for various values of p in the interval 0< p< 3 with 


-6 


an error within 5.1 x 10 ©. The introduction contains three chapters describing the construc- 
tion and use of the tables, and a fourth chapter, by Professor F. Mosteller, deals with appli- 
Cations. 
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